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I. 



INTRODUCTION 



A. BACKGROUND 

The United States Navy desires to provide standardized 
test facilities for the testing of its aviation gas turbine 
jet engines (Ref 1) . 

Today's high performance and high technology jet 
engines require accurate, consistent, and controlled 
testing to assess performance throughout the jet engine's 
entire operational range Modern day test facilities must 
incorporate state of the art equipment and instrumentation 
such that jet engines may be tested in or out of the 
aircraft, around the clock, and year round. 

B. NATURE OF THE PROBLEM 

Presently, 75 Navy and Marine Corp activities worldwide 
operate 81 Jet Engine Test Facilities (Ref 1) . Of these, 

60 are fully operational and 51 are over 25 years old. 
Immediate problem areas which hinder operational needs 
include: structural degradation, obsolescence, noise, and 

air pollution. 

1 . Structural Degradation 

High jet engine exhaust gas temperatures and water 
impingement on diffusers and interior concrete walls have 
caused cracking and spalling of test facility interiors. 
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Concrete reinforcing rods exposed to a hot, moist 
environment will corrode and thereby reduce structural 
integrity. Noise reduction baffles and exhaust gas turning 
vanes are being damaged by flying concrete debris. 
Replacement of baffling and turning vanes is required as 
often as every three years. The maintenance requirements 
navy wide has significantly increased the workload and 
operational costs of the Navy Public Works Department. 

2 . Obsolescence 

Several of the test facilities require major 
modernization or even complete replacement due to current 
aircraft designs and test requirements. Examples of 
modernization include improved engine handling systems and 
engine thrust stands. Increased access to engine 
components and test points greatly improves testing 
efficiency. Additionally, with today's high technology 
engines, state of the art data acquisition systems are 
required. 

3 . Noise 

Most of the older test facilities use water for 
exhaust gas cooling. The water and the high exhaust gas 
temperature deteriorates the acoutic baffling silencing 
capabilities. The increased . noise level degrades the work 
environment. Worker efficiency is reduced and the high 
test facility interior sound levels are being suspected of 
contributing to engine deterioration. State of the art 
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engines using thinner, lighter materials are more sensitive 
to acoustic fatigue (Ref 1) . 

4 . Air Pollution 

Today's high power and high mass flow engines 
constitute a major problem to the environment. Visible 
particulate matter is emitted in excess of allowable 
environmental pollution standards. By Executive Order 
Number 11282, May 26, 1966 (Ref 2), Jet Engine Test 
Facilities must meet local laws concerning the environment. 
Judicial action by the state of California against Jet 
Engine Test Facilities within California (Ref 3) has 
identified a need to address pollution control measures 
during the design, construction, and operation of Jet 
Engine Test Facilities. 

C. TYPES OF TEST FACILITIES 

The Navy operates five types of Jet Engine Test 
Facilities (Ref. 1) . Each facility utilizes some form of 
hold down device for the engine or the aircraft and 
instrumentation to monitor the performance of the engine. 
Each facility must be designed such that several engine 
types may be tested (e.g., turboprop, turbofan, turboshaft, 
turbojet) . The types of facilities are: 

1 . Run-up Pad 

The simplest test facility is the run-up (power 
check) pad. Normally, the engine is tested while installed 
in the aircraft. Therefore, the aircraft acts as the hold- 



down device. The run-up pad is usually located in remote 
locations due to high noise levels encountered during 
testing. 

2 . Sound Suppressors 

These are usually equipment packages which are used 
to reduce inlet and exhaust noise from a tied down 
aircraft. The sound suppressors may be used to augment or 
to substitute for the run-up pads. 

3 . Open Test Stand 

This is an outdoor test facility which supports and 
constrains the engine during test. Due to noise levels, 
the open test stand must be located in a remote area. 
Portable sensors and instrumentation are utilized to 
monitor engine performance. 

4 . Hush House 

This is a total aircraft enclosure. It permits 
round the clock and year round testing. Normally on twin 
engine aircraft, only one engine may be tested in full 
after burner. The other engine is held at idle. 

5 . Fully Enclosed Test Cell 

This test facility is normally capable of testing 
several types and sizes of jet engines. The engines may be 
in or out of the aircraft. The facility is normally 
equipped with all the necessary instrumentation and data 
gathering equipment required to perform routine or complex 
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testing after jet engine rework or overhaul. Testing may 
be performed round the clock and year round. 

The type of test facility to utilize at an activity is 
dependent upon several factors. Availability of military 
construction funds, required control of noise and emissions 
(environmental concern) , amount of maintenance testing, and 
military and civilian manpower requirements are all key 
factors . 

D. RESEARCH GOAL 

To provide quality jet engine testing, to comply with 
local environmental standards, and to ensure fleet 
readiness, the Navy embarked on a comprehensive program in 
1984 to upgrade and to modernize as necessary all jet 
engine test facilities under its cognizance. To this end, 
it was necessary to review existing jet engine test 
procedures and to examine jet engine test facilities to 
determine which facilities met desired criteria. 
Additionally, the Navy decided to explore the feasibility 
of the computer modeling of aviation gas turbine jet engine 
test facilities. It was felt that if a computer model 
could accurately represent an actual operational test 
facility then a computer model could also be used to 
predict the operation of future proposed test facilities. 
The purpose of this research was to computer model the 
Royal Danish Test Facility and to compare computer 
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generated temperatures, pressures, and velocities with 
actual data obtained in the field and to investigate 
techniques necessary to make computer simulations of flow 
and temperature fields useful for predictive purposes. 
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II. TEST FACILITY OPERATION 



The basic geometry of the Denmark test facility is 
shown in Figure C.l. For simplicity all storage areas, 
engine preparation rooms, and control modules are not 
included. The test facility is basically comprised of the 
primary air inlet, a test bay, an exhaust duct (also called 
an augmentor tube) , secondary air inlets, and an exhaust 
stack. Turning vanes, which are not shown, are located in 
the exhaust duct to exhaust stack transition region. This 
feature better utilizes the total exhaust stack cross 
section with an attendant reduction in exhaust gas noise 
(Ref 4) . Baffles (not shown) are located at three places: 
At the top of the exhaust stack, within the secondary air 
inlets, and within the primary air inlet (test bay doors) . 
The baffles also provide noise reduction (Ref 4) . 

During a jet engine test, the engine is aligned with 
the centerline of the exhaust duct. This can be 
accomplished either with the engine in the aircraft or with 
the engine supported in a test stand. With proper 
alignment, a near uniform flow profile develops (Ref 3) 
inside the augmentor tube. Primary air is drawn in through 
the primary air inlet (test bay door/baffle configuration) . 
For a F-100 Engine, at full throttle and after burner on, 
the primary air mass flow rate is 260 kg/sec of which 
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100 kg/sec is the mass flow rate through the engine 
(Ref. 4) . As the exhaust gas enters the augmentor tube, 
air from within the test bay is entrained into the 
augmentor tube. The exhaust gas moving through the 
augmentor tube also entrains air coming through the 
secondary air inlets. The secondary air mass flow rate is 
340 kg/sec (Ref. 4) . Both the entrained primary and 
secondary air play an important role in the test facility 
operation. Hot and noisy jet engine gases must be treated 
prior to leaving the test facility in order to comply with 
environmental standards. To treat jet engine exhaust, its 
temperature, noise level, and particulate density must be 
lowered . 

The ratio of total "cooling" air (primary and 
secondary) mass flow rate to the engine mass flow rate is 
called the Augmentation Ratio. The Augmentation Ratio is a 
key to successful test facility operation. With 
appropriate primary air and secondary air mass flow rates, 
the exhaust gas is sufficiently cooled which helps prolong 
duct/stack life and the exhaust gas is sufficiently diluted 
such that proper pollution levels are obtained (Ref. 2) . 

The Augmentation Ratio of the Denmark Test Facility is 5:1. 

With a computer model that can accurately predict test 
facility operation, design personnel will have the ability 
to easily change model parameters such as engine location 
or structural dimensions and to observe their effect on 
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test facility operation. For example by changing the 
secondary air inlet configuration, the augmentation ratio 
is changed. This in turn has an effect over exhaust gas 
cooling. Additionally, design engineers may wish to 
experiment with various exhaust stack heights or exhaust 
duct openings. A mere changing of a few lines of computer 
model source code accomplishes this. In all, a computer 
model could save significant time and money in the design 
of future jet engine test facilities. 
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III. DESCRIPTION OF PHOENICS 



A. WHAT IS PHOENICS 

PHOENICS is a general-purpose, computer-code system for 
simulating fluid flow, heat transfer, and chemical reaction 
processes in engineering equipment or the environment. 
PHOENICS is an acronym for Parabolic, Hyperbolic, or 
Elliptic Numerical Integration Code Series. Typical 
applications of PHOENICS include: Thermal hydraulics of 

nuclear reactors and other power plants; external 
hydrodynamics and aerodynamics of cars, ships, buildings, 
etc; air movement within an enclosed space, both under 
normal conditions and in the presence of a fire; flows in 
pumps and turbomachinery. PHOENICS was developed by D. 
Brian Spalding of CHAM Ltd and of Imperial College of 
London, England. Information regarding the history and 
development of PHOENICS is contained in the user's guide 
(Ref 5) . 

B. HOW PHOENICS WORKS 

The PHOENICS computer code consists of two separate 
computer programs: EARTH and Ql. The EARTH program 

resides in the computer mainframe as an object code and 
cannot be accessed by the user. The user couples his 
problem (defined by the user) to the EARTH code. The EARTH 



10 



program contains the main flow simulating software and 
various coding sequences which simulate the Laws of Physics 
as applied to the conservation principles: energy, mass, 

momentum. The Q1 program defines the problem and differs 
from one flow simulation program to the next. The Q1 
program represents the user's decision as to what problem 
is to be simulated. The Q1 program activates the EARTH 
program which solves the problem. 

The PHOENICS code is capable of describing flow 
simulations in terms of distributions in time and space. 

The distributions involve assigning numerical values for 
all dependent variables to an orderly array of cells 
located in the physical domain of the problem. 

Temperature, pressure, concentration, KE (Turbulence 
Kinetic Energy) , and EP (Kinetic Energy Dissipation Rate) 
values are located at a cell's center while velocity values 
are located at the positive faces of the cells. A cell 
configuration in the X-Y plane is illustrated in Figure 
C.2. Temperatures, pressures, and concentrations would be 
located at N, S, E, W, and G for example while velocities 
would be at N', S', E', and W' . A similar procedure is 
utilized for the Z (low-to-high) direction. This 
convention is called the "STAGGERED GRID" arrangement. If 
the flow simulation is time dependent, the distributions 
are obtained for each variable at a succession of times 
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instants. The Denmark test facility model flow simulation 
is for steady state conditions. 

In general, a PHOENICS flow simulation is four 
dimensional: 1 time and three space coordinates. For ease 

of reference to the user, PHOENICS measures time in the 
past-to-future direction; X coordinates in the west-to-east 
direction, Y coordinates in the south-to-north direction, 
and Z coordinates in the low-to-high direction. This 
concept is illustrated in Figure C.3. Although Figure C.3 
represents a Cartesian configuration, individual cells may 
vary in shape and size as determined by a user defined 
code. This option, called BFC, was not utilized in this 
investigation. The Denmark test facility model is a 
Cartesian configuration of cells which vary in size. 

PHOENICS calculates the dependant variables over the 
complete set of cells, which defines the operational (e.g., 
system, physical) domain, with linear algebraic equations 
that represent the conservation laws. Here the partial 
differential equations are those governing the conservation 
of energy, momentum, and mass and the linear algebraic 
equations are obtained by the finite volume method. The 
finite volume method is a method having some features from 
the finite differencing method and the finite element 
method. The conservation partial differential equations 
will not be discussed here; except to say that the original 
non linear, coupled conservation partial differential 



equations are linearized. Interested readers are directed 
to the PHOENICS User's Manual (Ref. 5) for detailed 
information. 

PHOENICS utilizes an iterative process involving a 
multi-stage sequence of dependent variable value 
adjustments. The process is repeated several times until 
the left and right imbalance (called a residual) of each 
equation is reduced to a negligible amount. A variety of 
iterative procedures are available to the user. The 
"SLAB-SWEEP" iterative process feature of PHOENICS for 
equation imbalance reduction was chosen for the Denmark 
test facility model. Slabs are arrays of cells having the 
same Z (low-to-high) coordinate. PHOENICS conducts many 
mathematical operations over a slab in an attempt to reduce 
dependent variable imbalances among cells within the slab 
prior to moving on to the next slab. The adjustment 
operations are referred to as "SLABWISE SOLUTIONS". The 
operations begin at the lowest numbered slab and continue 
sequentially to the highest numbered slab. A sweep is a 
completed set of slabwise operations. The equations for 
dependent variable values at one slab are coupled to values 
of other dependent variables at neighbor slabs. For this 
reason numerous sweeps are normally required until all 
dependent variable imbalances are negligible. This is a 
necessary but not a sufficient condition for problem 



convergence. 



C. PHOENICS INPUT LANGUAGE 

The source code used in a PHOENICS Q1 program is 
arranged in a 24 group data structure. Figure C.4 lists 
the group categories of the data structure. The developers 
of the PHOENICS code suggest this ordering of the groups 
normally conforms to a well thought out flow simulation 
problem description. However, any group ordering the user 
desires will be acceptable to PHOENICS. It is the PHOENICS 
defined variables and keywords plus the user defined 
variables that actually trigger desired PHOENICS responses 
as opposed to a particular group ordering. A detailed 
description of group functions is explained in the PHOENICS 
User's Manual (Ref. 5) . 

D. MODEL GEOMETRY 

PHOENICS utilizes three distinct procedures for 
defining a flow simulation operational domain (grid) : 
Cylindrical-polar, Body-fitted, and Cartesian coordinate 
systems . 

1. Cylindrical-Polar Coordinate System 

The cylindrical-polar grid arrangement consists of 
planes with constant Z coordinates that are perpendicular 
to an axis of rotation, planes with a constant X coordinate 
(usually called the Theta coordinate) that intersect with 
the axis of rotation giving the X coordinate a 
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representation in radians, and concentric cylindrical 
surfaces with a constant Y radial coordinate. 

2 . Body-Fitted Coordinate System 

This arrangement is best thought of as starting 
with an orthogonal group of cells made of clay. The cells 
are then squeezed, bent, twisted, and stretched into the 
desired shape for the fluid flow simulation operational 
domain. Each cell although distorted in shape and size 
still remains in contact with its original neighbors. 

3 . Cartesian Coordinate System 

The Cartesian grid arrangement is composed of cells 
in a mutually perpendicular X, Y and Z coordinate system. 

As such the cell faces form planes that intersect at right 
angles. The spacing between the cell faces (planes) is 
user defined. 

Due to the regularity of a Cartesian coordinate 
system, the information required to specify the geometry is 
rather modest. The PHOENICS variables NX, NY, and NZ 
specify the number of cells the operational domain is to 
have in the X, Y, and Z directions respectively. The 
PHOENICS variables XULAST, YVLAST, and ZWLAST specify the 
domain lengths in the X, Y, and Z directions. The "METHOD- 
OF-PAIRS" is a PHOENICS coding sequence to specify the 
number of cells per length of domain in the X, Y, and Z 
directions. XFRAC, YFRAC, and ZFRAC are PHOENICS variables 
for arrays used in the method-of-pairs coding sequence. 



The odd indexed array elements, for example, XFRAC(l), 
XFRAC(3) , XFRAC(5), etc., specify the number of cells in a 
particular length of operational domain while the even 
indexed array elements XFRAC(2), XFRAC(4), and XFRAC(6) 
specifies the length of each cell. A negative argument for 
the first element in each array tells PHOENICS that the 
method-of-pairs is being utilized. For example, "XFRAC(l) 
=-2 . 0 ; XFRAC (2 ) =0 . 1" tells PHOENICS the method-of-pairs is 
to be used and that the first cell group in the X direction 
is to have two cells and the length of each cell is 0.1 
meter. "XFRAC (9) =4 . 0 ; XFRAC (10) =0 . 15" specifies the number 
of cells and lengths of each cell in the "FIFTH" X 
direction cell group. 

E. POROSITY AND BLOCKAGE 

In any flow simulation problem there exists walls, 
submerged bodies, or some form of obstruction to fluid 
flow. For example, in Hagen-Poiseuille pipe flow the fluid 
is confined to the interior of the pipe. Additionally, due 
to wall friction the fluid's velocity at the wall is zero. 
In a PHOENICS flow simulation problem, obstructions to 
fluid flow must also be introduced. In a jet engine test 
facility model one must be concerned with obstructions to 
flow as a result of walls, floors, baffle material, turning 
vanes, and the jet engine itself. To account for the 
physical obstructions, PHOENICS makes use of the command 
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"CONPOR (R. R, TYPE , XF, XL, YF ,YL,ZF,ZL)". The real number 
"R.R" normally has a value from 0.0 to 1.0. This value 
tells PHOENICS how much obstruction to flow is to be 
introduced by a surface or volume. For maximum obstruction 
(blocked flow), "R.R" would be 0.0. The "TYPE" argument 
specifies a surface or volume. For example a surface might 
be the inside of the test bay ceiling. A volume might be 
the geometric body of the jet engine. For a surface 
obstruction the "TYPE" Arguments "WEST, EAST, SOUTH, NORTH, 
LOW, OR HIGH" (the faces of an individual cell) are used. 
For a volume obstruction the type "CELL" is used. The 
arguments "XF, XL, YF, YL, ZF, and ZL" identify the first 
and last cells of the obstruction domain in the X, Y, and Z 
directions respectively. For example "CONPOR (0.0, HIGH, 
1,3, 1,1, 3, 3)" tells PHOENICS to introduce complete 
obstruction to flow through the "HIGH" (Z direction) 
surfaces of the group of cells starting at IX=1, IY=1, and 
IZ=3 , and ending at cells IX=3, IY=1, and IZ=3 
respectively. Each cell has an IX, IY, and IZ identity. 
Figure C.5 illustrates "HIGH" surface obstruction to flow. 
The cross-hatched surface blocks Z direction flow. "CONPOR 
(0. 0, CELL, 1, 3 , 1, 1, 1, 3) " would introduce maximum obstruction 
to flow through the surfaces of the volume defined by these 
coordinates. In Figure C.5, the "CELL" volume outlined 
with bold lines prevents fluid flow perpendicular to any of 



its surfaces. 



PHOENICS also has the ability to introduce wall or 
surface "FRICTION" which brings about zero magnitude 
velocity components on a cell's surface and calculates the 
friction coefficient. For the above "HIGH" example, the 
magnitude of flow in the Z direction was set to zero on the 
surfaces of several cells. To also set the X and Y 
velocity components to zero on the same "HIGH" surfaces, 
the CONPOR command requires a slight modification. 

Utilizing a "NEGATIVE" argument with the CONPOR command 
accomplishes the desired action. For example 
"CONPOR (0.0, HIGH, 1,3, 1,1, 3, -3)" causes the X, Y, and Z 
velocity components on the "HIGH" surfaces of that 
particular group of cells to be zeroed. Through the use of 
the negative argument convention, velocity components on 
the external surfaces of a volume (cell) can be zeroed, and 
friction introduced into the model. 

F. SOURCE AND BOUNDARY CONDITIONS 

In any PHOENICS flow simulation problem description, 
the user's Q1 program must introduce source and boundary 
conditions to the EARTH program. Through the use of 
PHOENICS keywords, source and boundary conditions may be 
identified as energy, momentum, or mass flow. Energy 
sources are identified with the enthalpy variable "H" , 
momentum by the velocity variables "U, V, or W" , and mass 
flow utilizes the pressure variable "P". By using the 
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notation "HI” a fluid's first phase enthalpy would be 
identified while "H2" would refer to second phase enthalpy 
for a two phase flow simulation. Information about inflow 
and outflow of any variable in the system domain is 
conveyed to the EARTH Program by way of the PHOENICS 
combination "PATCH-COVAL" statements. The PATCH command is 
used to specify the location at which a source or boundary 
condition is to be applied and to also specify the "TYPE" 
factor or what the PATCH command is to do to the system 
domain. The COVAL command is used to specify the 
coefficient and value which combine to determine the 
magnitude of the dependent variable over a PATCH. For 
example, "PATCH (JETOUT, LOW, 1, 1, 5 , 6 , 8 , 8 , 1, 1) ; 

COVAL (JETOUT, PI, FIXFLU, 100. 0) ; 

COVAL (JETOUT , W1 , ONLYMS ,500.0); COVAL (JETOUT , HI , ONLYMS , 
200000.0)" specifies a 100 kg/sec, steady-state, fixed-flux 
mass source. The mass source, named JETOUT, begins at the 
low faces of the cells IX=1, IY=5, and IZ=8 and ends at 
IX=1, IY=6, and IZ=8. The mass source has a Z direction 
velocity of 500 m/sec and an enthalpy of 200,000 J/kg. The 
COVAL coefficient "ONLYMS" specifies that only convection 
flow effects are influential in the transport of enthalpy 
from the cell faces to the cell centers. That is, 
diffusive effects are zero. Appendix 2 of the PHOENICS 
User's Manual (Ref. 5) lists the various "TYPES", 
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"COEFFICIENTS", and "VALUES" that may be used with the 
combination "PATCH-COVAL" statements. 

G. INITIAL CONDITIONS 

In order to obtain the solution of a steady state 
problem, PHOENICS starts the iterative solution with an 
initial field of the dependent variables. The default 
initial fields for all dependent variables is zero over the 
entire physical domain. However, convergence to the true 
solution can be speeded up by the user specifying initial 
field conditions via "FIINIT" (Field Initialization) 
statements over the entire operational domain or 
"PATCH-INIT" statements over subdomains of the system. For 
example "FIINIT (HI) =AMOUNT" initialize a first phase 
enthalpy field to a user-defined value "AMOUNT". 

"PATCH (GAS IN , INIVAL, XF, XL, YF , YL, ZF , ZL, TF , TL) ? INIT (GASIN, PI , 
COEFFICIENT, VALUE)" initializes a first phase pressure 
field for the user-defined PATCH "GASIN". XF, XL, YF, YL, ZF, 
and ZL define the geometric region of the PATCH called 
"GASIN". 

H . MAXIMUM/MINIMUM VALUES 

At the beginning of a flow simulation problem iterative 
process, a poorly defined dependent variable field value 
may have an adverse affect on problem convergence. At 
times, dependent variables "STRAY" outside of a physically 
meaningful range and cause numerical divergence. To help 
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minimize divergence of the solution PHOENICS utilizes the 
"VARMIN" and VARMAX" commands. These commands are used to 
set upper and lower bounds on any dependent variables. 
"VARMIN" (VI) =SLOW; VARMAX (W1)=FAST" sets a lower limit of 
"SLOW" to VI and an upper limit of "FAST" to Wl. 

I . UNDER-RELAXATION FACTORS 

A feature of PHOENICS which greatly assists problem 
convergence is the application of "UNDER-RELAXATION" 
factors to flow simulation dependent variables. This 
practice is most useful on an operational domain comprised 
of ill-proportioned cells (very large cells adjacent to 
very small cells) . During program execution, EARTH'S 
dependent variable correction equations may develop 
enormously large coefficients between adjacent cells. The 
coefficient disparity often contributes to solution 
divergence. With selective "UNDER-RELAXATION" applied to a 
flow simulation program, PHOENICS introducers " FALSE-TIME - 
STEP" values to select dependant variables thereby 
preventing large dependent variable changes in magnitude 
per sweep. The "FALSE-TIME-STEP" brings equation 
coefficients closer together "LESS ABRUPTLY" which promotes 
problem convergence. "UNDER-RELAXATION" for pressure, 
velocity, and enthalpy components is frequently required 
for steady state problems. 
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J. PROGRAM RUN PROCEDURES 

The Denmark test facility model was analyzed on the 
Naval Postgraduate School IBM-3033 mainframe computer 
system. Program run procedures involve interactive 
terminal sessions and a batch submission. The necessary 
flow simulation Q1 program and Job Control Language (JCL) 
programs reside in the User's mainframe storage area. All 
flow simulation programs are submitted for execution by 
means of "EXEC" (execute) and JCL files. The EXEC files 
are activated by the keywords "C0PYQ1", "RUNSAT", and 
"RUNEAR" . 

The C0PYQ1 EXEC file queries the user as to the name 
and file type of the flow simulation Q1 program source 
code. Afterwards the C0PYQ1 EXEC file stores the source 
code and appropriate JCL in a temporary file. At this 
point, the user must ensure the above transfer of data was 
successful. This is accomplished through the user's 
"VIRTUAL READER" file. An inspection of this file will 
reveal various condition codes which indicate the status of 
data transfer. Figure C.6 is an example of C0PYQ1 EXEC and 
JCL files. 

The "RUNSAT" EXEC file extracts from the user's 
mainframe storage area JCL which enables a Q1 program to 
EARTH program communications link. This JCL is also stored 
in a temporary file. Again a virtual reader file 
inspection must be made to ensure that a successful data 
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transfer was achieved. Figures C.7 and C.8 are examples of 
RUNSAT EXEC and JCL files. 

Finally after successful C0PYQ1 and RUNSAT data 
transfers, the flow simulation problem is ready for 
execution. The "RUNEAR" EXEC file is called, Necessary 
mainframe storage is allocated and the PHOENICS system is 
activated. Figure C.9 is an example of RUNEAR EXEC and JCL 
files. Upon program completion, PHOENICS output is sent to 
the user's virtual reader file. After the PHOENICS output 
is analyzed, the user may discard the information or 
utilize a PHOENICS restart which allows the problem to 
continue from where it last ended. Upon completion of a 
PHOENICS run, output is stored on disk file "DF09" . To 
access this file and to continue the problem the PHOENICS 
code "RESTRT (A, B, . . . . , Z) " must be utilized in the Q1 
program. "A,B,...,Z" represent the dependent variables 
which are to continue from their last value when the 
program ended. If the user desires all dependent variables 
to restart from their last values then "RESTRT (ALL) " can be 
used in the Q1 program. 
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IV. MODELING THE DENMARK TEST FACILITY 



A. MODEL DESCRIPTION 

The Denmark test facility model utilizes a Cartesian 
grid arrangement. Figure C.10 is a schematic diagram of 
the test facility in the Y-Z pane (side view). Figure C.ll 
is the X-Z plane (top view) representation. The PHOENICS 
code allows for a plane of symmetry. For the Denmark test 
facility model the vertical Y-Z plane is one of symmetry. 
Therefore, in Figure C.ll only half of the Denmark test 
facility is represented. In this figure the jet engine is 
slightly of f-of-center from the Z axis for clarity only. 

To facilitate the Denmark test facility model geometry 
specifications, the system domain was divided into user 
defined zones. Each zone comprises a particular section of 
the Denmark test facility model. For example in Figure 
C.12, "YLBL1" is a user defined variable which specified 
the length of the first zone in the Y direction. This 
particular zone represents the length in the Y direction 
from the ground or floor of the Denmark test facility up 
2.0 meters. n ZLEXT2" represents the length in the Z 
direction of the atmosphere behind the test facility, and 
so on. Figure C.12 is a schematic diagram of the system 
domain zones which make up the vertical Y-Z plane while 
C.13 illustrates the zones in the horizontal X-Z plane. 
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The zone concept has proven to be very flexible. The 
number of cells in each zone may be specified by a user 
defined variable. For example, "XNCBL1" specifies the 
number of cells in the X direction for the length of 
operational domain "XLBL1". "ZNCEX2" specifies the number 
of cells in the Z direction for the length of operational 
domain "ZLEXT2" and so on. Grid refinement for a 
particular zone can be accomplished without any difficulty 
by changing one variable in the Q1 code. Figure C.14 is 
the method-of-pairs coding sequence which specifies the 
operational domain zones and numbers of cells used in the 
Denmark test facility model. Figure C.15 is an example of 
the EARTH output for the Q1 code in Appendix A. 

Program run time, program cost and EARTH output are 
directly affected by the number of cells in the model of 
the physical domain. By increasing the number of cells in 
the X, Y, or Z direction zones, the EARTH program requires 
more "SLAB-WISE" Computations thereby increasing computer 
CPU requirements and costs. By increasing the number of 
cells in the X or Y directions (i.e., NX or NY), the number 
of cells and therefore the number of calculations is 
increased. For example for a fixed NY, SAY 16, an increase 
of 1 cell in the X direction will result in 7 (number of 
dependent variables) times 16 = 102 additional equations 
(and unknowns) . An increase of an additional cell in the Z 
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direction will result in 7 times NX times NY additional 
equations. For example for NX=9 and NY=16, increasing NZ 
from 28 to 29 will result in an additional 1,008 equations. 
Therefore, deciding to use a "COURSE" or "FINE" model of 
the physical domain is determined by the user's time and 
financial constraints and the desired information from the 
flow simulation program. The Denmark test facility model 
has NX=9, NY=16 and NZ=28 for a grand total of 4,032 cells. 
Since each cell has 7 dependent variables, 28,224 linear 
algebraic equations must be solved iteratively to achieve 
problem convergence. This is a big job even for an IBM 
mainframe computer. 

For the purpose of this research, the Denmark test 
facility was modeled with a "COARSE" grid with the intent 
of achieving problem convergence prior to using a "fine" 
grid. Due to time constaints only one grid was used in 
this research. Therefore, problem convergence could not be 
identified as "GRID INDEPENDENT". 

Although the Denmark test facility model grid is 
considered "COARSE", many of the cells vary in size. The 
number and size of cells for a particular region of the 
grid was chosen based upon the expected activity in a 
region. For example, at the primary air inlet, one would 
expect primary cooling air to be entering at a relatively 
low velocity and in a laminar and uniform flow. Then as 
the primary air approaches the jet engine intake, fluid 
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velocities should be increasing in magnitude and changing 
in direction. At the outlet of the jet engine there exists 
a very hot and a very fast fluid flow mixing with entrained 
primary cooling air. In this particular region with a lot 
or activity, a less "COARSE" grid is more desirable in 
order to incorporate the mixing process into the model. As 
the exhaust gas enters the augmentor tube it mixes with the 
cooler entrained secondary air. Again, a less "COARSE" 
grid is necessary to better analyze the activity involved. 

Figure C.16 is a schematic diagram of the system domain 
grid for the X-Z plane while Figure C.17 represents the Y-Z 
plane. The two figures illustrate a less "COARSE" grid 
arrangement at the jet engine exhaust, augmentor tube 
inlet, and secondary air inlet. Again, grid refinement for 
any physical domain zone is accomplished by changing one 
user defined variable in the Q1 code. 

Although the Denmark test facility was modeled as a 
three-dimensional structure. The body fitted coordinate 
option of PHOENICS was not utilized. As such, the 
cylindrical jet engine, the augmentor tube, and the exhaust 
stack turning vanes were modeled as rectangular cubes. The 
jet engine was specified as a solid block with a hot mass 
source at one end and a mass sink at the other end. The 
jet engine was modeled with CONPOR and PATCH-COVAL 
statements . 
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The Denmark test facility walls, floor, and ceilings 
were modeled utilizing CONPOR and PATCH-COVAL statements. 
The CONPOR statements were given a zero porosity value 
thereby preventing mass transfer across any solid 
structure. PATCH-COVAL statements were necessary for the 
walls, floor, and ceilings in order to introduce 
appropriate friction and set fluid velocity components to 
zero magnitude at the structure surfaces. 

B. MODEL BOUNDARY CONDITIONS AND TURBULENCE MODELS 

Boundary conditions for the Denmark test facility 
openings to the atmosphere were specified with PATCH-COVAL 
statements. Ambient pressure was set to 101,325 Pascals 
while ambient temperature was set to 295 K. 

To incorporate turbulence into the flow simulation 
problem, this research utilized the PHOENICS K-E turbulence 
model. The K-E model introduces the KE (turbulence kinetic 
energy) and EP (kinetic energy dissipation rate) dependent 
variables into the problem. Turbulence modeling is not an 
exact science and the choice of the K-E model for this 
research was arbitrary. The PHOENICS user's manual 
(Ref. 5) explains in detail the K-E turbulence model and 
three other turbulence models available to the PHOENICS 
user. 

Readers may refer to Appendix A, the annotated Q1 code, 
which explains in detail the user defined variables for the 
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Denmark test facility model. Columns 1 and 2 are reserved 
for the beginning of executable statements, while non- 
executable statements (comments) begin in column 3 or 
beyond. 
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V. ANALYSIS OF RESULTS 



A. OBTAINING A FINAL SOLUTION 

The original purpose of this research was to computer 
model the Royal Danish Air Force Gas Turbine Test Facility 
using the PHOENICS code. The computer model would then be 
used to compare field data obtained from the test facility 
with computer generated pressures, temperatures, and 
velocities. If the computer generated output compared well 
with the field data, one might suggest that the computer 
model was successful. However, if field data from a test 
facility was non-existent or perhaps very limited how then 
could the success of a computer model be measured? 

The field data from the Denmark test facility was 
limited in scope. Besides the mass flow described in 
Chapter II, the only other field data available are the 
test bay pressure and the velocity and temperature of the 
engine exhaust at one point in the augmentor tube and at 
one point in the exhaust stack. Since the field data from 
the Denmark test facility was limited, a comprehensive 
comparison of field data with numerical results was not 
possible. To investigate the EARTH output (numerical 
results) and conclude that the numerical results were 
representative of the test facility, this study examined 
the EARTH output from three viewpoints: common sense 
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approach, conservation equation balances, and problem 
convergence. 

B. COMMON SENSE APPROACH 

The Denmark test facility was modeled with five major 
regions or zones. Within each zone a common sense approach 
to what one would expect the gas pressure, temperature, and 
velocity conditions to be was taken. The five major model 
zones are: primary air inlet, test bay, secondary air 

inlet, augmentor tube, and exhaust stack. The EARTH output 
presented in Appendix B is the numerical results of only 5 
slabs in the major model zones. 

1 . Primary Air Inlet 

The primary air inlet is comprised of the following 
cells: IX=1 to IX=8 , IY=1 to IY=12, and IZ=2. Cool ambient 
air (295 K) is drawn into the test facility through the 
primary air inlet (test facility doors) . Pressure drops at 
the test facility doors were very small with a maximum drop 
of about 12 Pascals. Air flow through the primary air 
inlet was predominantly in the Z-direction (W) . Velocities 
ranged from 1 to 2 m/sec with an average of about 1.7 
m/sec. The horizontal X-direction velocity (U) and the 
vertical Y-direction velocity (V) components were 
negligible. 

With an average gas density of about 1.2 kg/m at 
the test facility doors, a primary air mass flow rate of 
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about 135 kg/sec was calculated. Recall from chapter IV 
that the Q1 program modeled only one half of the test 
facility due to the Y-Z plane of symmetry. Therefore, the 
135 kg/sec primary air mass flow rate is for one half of 
the structure. The calculated primary air mass flow rate 
for the entire structure would be about 270 kg/sec. Of 
this amount, about 170 kg/sec would be primary cooling air 
and the remaining 100 kg/sec flows through the engine. 

From this point on, references to mass flow rates will be 
for the entire test facility. The temperature at the test 
facility doors remained at about 295 K. At the test 
facility doors, one would not expect much activity and a 
laminar flow. The KE and EP turbulence dependant variables 
were negligible which indicated laminar flow. 

2 . Test Bay 

The test bay comprises cells IX=1 to IX=8, IY=1 to 
IY=12, and IZ=2 to IZ=9. The jet engine which was modeled 
as a 2000 K, 100 kg/sec sink-source is located at cells 
IX=1, IY=5 to IY=6 , and IZ=6 to IZ=7. At 2000 K, the 
density of the engine exhaust gas was 0.18 kg/m 3 . The jet 
engine provides the excitation which draws ambient air in 
through the primary air inlet, entrains primary cooling air 
into the augmentor tube from the test bay and entrains 
secondary cooling air into the augmentor tube from the 
secondary air inlets. As primary air moved through the 
test bay from the primary air inlet to the sink side of the 
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jet engine, its pressure and velocity changed. At the two 
cells (IX=1, IY=5 to IY=6, and IZ=5) before the engine 
inlet, about 1.5 meters upstream of the inlet, primary air 
pressure dropped about 40 Pascals. The primary air W- 
velocity component was about 3 m/sec at IZ=4, about 3 
meters upstream of the engine inlet. The horizontal U- 
velocity component was about -4.5 m/sec which indicated 
that air was being drawn into the engine from the side 
walls of the test bay. The vertical V-velocity component 
was about -4 m/sec above the engine and about 4 m/sec below 
the engine. This indicated that air was being drawn down 
from the test bay ceiling into the engine and at the same 
time being drawn up from the test bay floor. The turbulent 
parameters KE and EP were very small indicating laminar 
behavior at the engine inlet. The temperature at the 
engine inlet remained at about 295 K. 

At the discharge side of the engine (cells IX=1, 

IY=5 to IY=6, and IZ=8) , significant activity was observed. 
The pressure increased about 2100 Pascals. In adjacent 
cells above/below and longside the engine at IZ=8, the 
pressure dropped about 410 Pascals. The pressure 
difference among the cells caused an expanding or 
diverging effect if the jet engine exhaust gas. The U- 
velocity component of the gas was about 30 m/sec while the 
V-velocity component was about 50 m/sec. The jet engine 

exhaust gas expanded out towards the rear wall of the test 
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facility as opposed to going directly into the augmentor 
tube. However, due to the large magnitude of the W- 
velocity component, about 430 m/sec, the resultant 
direction of the engine exhaust was primarily into the 
augmentor tube. Hot, fast moving gas at the engine outlet 
mixed with the cool, slow moving entrained primary cooling 
air. Due to this mixing, small turbulence resulted. this 
premise is supported by the magnitude of the KE and EP 
parameters in the region. KE had a maximum value of about 
110 J/kg, while EP had a maximum value of about 5400 
J/kg/sec. Due to this mixing of engine exhaust and primary 
cooling air, over a distance of about 1.5 meters, the 
exhaust gas was cooled from 2000 K to about 1830 K prior to 
entering the augmentor tube (IZ=10) . Over the same 1.5 
meter distance, the exhaust gas W-velocity component 
decreased from about 430 m/sec to about 400 m/sec. 

The investigation of the test bay was localized to 
a small region around the jet engine. For the most part 
the activity in the rest of the test bay was uneventful 
(i.e. no turbulence and insignificant pressure, 
temperature, or velocity changes) . 

The above observations in the test bay and at the 
inlet and outlet ends of the jet engine are realistic. One 
would expect primary air velocity components and pressure 
to change somewhat prior to going into the engine and more 
significantly at the engine outlet. Additionally, due to 
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the primary air being drawn into the engine a low pressure 
area at the inlet makes good sense. While at the outlet 
side of the engine a high pressure is expected. At the 
engine outlet one would expect that the hot, fast engine 
exhaust colliding and mixing with the cool, slow entrained 
primary cooling air would result in some turbulence. 

3 . Secondary Air Inlet 

As was mentioned earlier in chapter II, the 
secondary cooling air plays a significant role in the 
operation of a jet engine test facility. The secondary 
cooling air serves the purpose of providing additional cool 
(ambient) air for mixing with the hot engine exhaust. For 
the Denmark test facility the largest mass flow rate occurs 
through the secondary air inlets. The test facility data 
indicated a mass flow rate of about 340 kg/sec. In the 
structure model, secondary cooling air enters the augmentor 
tube at cells IX=5, IY=1 to IY=10, and IZ=10 to IZ=12. The 
secondary cooling air had an average U-velocity component 
of about -14 m/sec, which indicated that secondary cooling 
air was being drawn into the augmentor tube from outside 
the structure. 

A maximum pressure drop of about 480 Pascals at the 
secondary air inlet was in consonance with the relatively 
high U-velocity component. With an average air density of 
about 1.0 kg/m 3 at the secondary air inlet, a mass flow 
rate of about 420 kg/sec was calculated. Although the 
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secondary air inlet had a significant pressure drop and a 
relatively large U-velocity component (about -14 m/sec) it 
did not demonstrate turbulent behavior based upon the 
relatively low KE and EP parameters. The temperature at 
the secondary air inlet indicated a temperature rise of 
about 15 K. The results showed that warm air from the 
exhaust stack was flowing into the secondary air inlet 
region . 

With a primary cooling air mass flow rate of about 
170 kg/sec and a secondary cooling air mass flow rate of 
about 420 kg/sec, an Augmentation Ratio of about 5.9:1 was 
calculated . 

4 . Augment or Tube 

Three regions of the augmentor tube were chosen for 
investigation: The region IZ=10 to IZ=16, the slab IZ=20, 

and the slab IZ=22. At the beginning of the augmentor 
tube, IZ=10 to IZ=16, a large mass of primary and secondary 
cooling air collided and mixed with the hot, fast exhaust 
gas. As a result, significant turbulence occurred. A 
maximum KE value of about 7,100 J/kg was observed while 
maximum EP was about 500,000 J/kg/sec. In the 8.5 meters 
that the engine exhaust traveled, it cooled from 2000 K to 
about 920 K. The engine exhaust W-velocity component 
decreased from about 430 m/sec to about 125 m/sec. The U 
and V-velocity components were considered negligible. 
Figures C.18 and C.19 are PHOENICS enthalpy (HI) and 
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velocity (Wl) profiles for the cells IX=1, IY=1 to IY=10, 
and IZ=16. The reader is reminded at this point that an 
enthalpy value divided by a specific heat value of about 
1000 J/(kg K) will give a good temperature approximation 
within a cell. The profiles indicated minimum temperature 
and W-velocity values at the bottom (IY=1) and top (IY=10) 
of the augmentor tube and maximum temperature and W- 
velocity values at the centerline of the augmentor tube. 
These temperature and W-velocity distributions make good 
sense. Since the location IY=5 is at the height of the 
engine outlet, one would expect the engine exhaust to be 
hotter and to be moving faster here than at the top/bottom 
of the augmentor tube. Figures C.20 and C.21 are 
temperature and W-velocity profiles for the cells at IX=1 
to IX=5, IY=5 , and IZ=16. Although these profiles are at 
the same slab location (IZ=16) as figures C.18 and C.19, 
the examination of the variables is from the center of the 
augmentor tube to the side wall. These distributions 
indicated maximum values of the variables occurred at the 
center of the augmentor tube and minimum values at the side 
walls. These profiles also make good sense. In the cells 
at the augmentor wall (IX=5, IY=1 to IY=10, and IZ=16) , the 
secondary cooling air entering the augmentor tube has only 
mixed with the engine exhaust for about a 3 meter distance 
down the augmentor tube. 
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At augmentor tube location IZ=20 (about 22 meters 
downstream of the engine outlet) , the total mass flow was 
examined. The W-velocity component ranged from about 25 
m/sec to 75 m/sec with an average of about 45 m/sec. The U 

and V-velocity components were considered to be negligible. 

• • ^ 

With an average gas density of about 0.6 kg/m , the total 

mass flow rate was about 700 kg/sec. At augmentor tube 

location IZ=20, there was still some turbulence, which 

indicated continued mixing of the primary/secondary cooling 

air and the engine exhaust. The maximum KE value was about 

740 J/kg while EP maximum was about 10,400 J/kg/sec. The 

temperature ranged from about 520 K to about 635 K with an 

average of about 550 K. Figures C.22 and C.23 are 

temperature and W-velocity profiles for the cells at IX=1, 

IY=1 to IY=10, and IZ=20. Figures C.24 and C.25 are 

temperature and W-velocity profiles for the cells at IX=1 

to IX=5, IY=5, and IZ=20. A more uniform distribution of 

temperature and W-velocity was observed across the slab at 

IZ=20 than at slab IZ=16. This is reasonable since the 

total mass of primary and secondary cooling air and engine 

exhaust has mixed for a distance of about 22 meters 

downstream of the engine outlet. 

In the augmentor tube/exhaust stack transition 
region (IZ=22) , the test facility field data indicated the 
W-velocity component on a centerline with the jet engine 
was about 60 m/sec and the average temperature was about 
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650 to 675 K. The numerical results indicated a W-velocity 
component of about 63 m/sec and a temperature of about 
615 K. In a distance of about 4 meters from IZ=20 to 
IZ=22, turbulence decreased by about 50 percent. KE was 
about 400 J/kg while EP was about 4,500 J/kg/sec. 

5 . Exhaust Stack 

Field data from the Denmark test facility indicated 
a temperature rise of about 20 K from the bottom of the 
exhaust stack to the area under the exhaust stack baffle 
material. In the numerical model a temperature drop of 
about 50 K was observed. The Denmark test facility was 
modeled without baffle material at the top of the stack. 

Gas velocities at the top of the computer model stack were 
for the most part in the V-direction. The V-velocity 
ranged from about 10 to 45 m/sec with an average of about 
20 m/sec. Turbulence at the top of the exhaust stack was 
considered to be negligible. Average temperature at the 
top of the exhaust stack was about 550 K. Figures C.26, 

C . 27 , C. 28 , and C . 29 are temperature and V-velocity 
profiles for the top (IY=13) of the exhaust stack. Figures 
C.26 and C.27 are for the cells IX=1, IY=13, and IZ=21 to 
IZ=26 (along the Z axis). Figure C.26 indicated highest 
temperature values towards the rear (IZ=25) of the exhaust 
stack. Additionally, Figure C.27 indicated that the 
V-velocity fluctuated (low to high to low ...) from IZ=21 
to IZ=26. It is suggested that the manner in which the 
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turning vanes were computer modeled caused a lack of 
uniformity in the temperature and V-velocity distribution 
for this stack region. 

Figures C.28 and C.29 are the temperature and V- 
velocity distributions for the cells at IX=1 to IX=5, 

IY=13, and IZ=23 (from the center of the exhaust stack out 
to the side wall) . These figures indicated fairly uniform 
temperature and V-velocity profiles throughout the cells in 
this region. 

C. COMPUTER DATA AND FIELD DATA COMPARISON 

The difference between field data results and numerical 
results can be attributed to several things: baffle 

material (not included in the numerical model) , turning 
vanes (crudely modeled as rectangular cubes) , turbulence 
model (modeling is not an exact science) , and model grid 
(grid independence not verified) . 

The computer model total mass flow rate of about 700 
kg/sec is about 15 percent larger than the 600 kg/sec total 
mass flow rate of the Denmark test facility. It is 
suggested that not including the baffle material, in the 
computer model primary and secondary air inlets and the 
exhaust stack, resulted in larger numerical mass flow rates 
than the test facility mass flow rates. Test facility 
baffle information was not available for this research. 



Field data indicated a maximum pressure drop in the 
test bay of about 500 Pascals. The computer model 
calculated a maximum pressure drop of about 410 Pascals in 
the cells adjacent to the jet engine outlet. 

As was mentioned earlier, the exhaust gas temperature 
and W-velocity component in the augmentor tube/exhaust 
stack transition region compared very well with the field 
data . 

D. CONSERVATION EQUATION BALANCES 

The PHOENICS output code enables the investigator to 
determine the degree conservation of mass has been 
achieved. The PHOENICS output symbol for a mass flow 
source is "Rl" . By an algebraic summing of the mass flow 
sources, a conservation of mass flow calculation can be 
made. For a total accounting of all the mass flow sources 
in a computer model, the algebraic sum of the sources 
should be zero. For the Denmark test facility this sum was 
equal to 0.02 kg/sec. If this amount is compared to the 
100 kg/sec mass flow through the jet engine there is a 0.02 
percent error. Compared to the 700 kg/sec mass flow 
calculated for the augmentor tube the percent of error is 
reduced to 0.003. Based on this analysis, a mass balance 
or a continuity of flow exists in the Denmark test facility 



model. A mass balance or a continuity of flow is a 
necessary but not sufficient condition for problem 
convergence . 

E . CONVERGENCE 

Chapter III described the iterative process which 
adjusts the values of the dependent variables during the 
solution process. During the process, imbalances occurred 
between the left and right hand side of the dependent 
variable discretized algebraic equations. These imbalances 
are called residuals. Problem convergence would be 
obtained when all residuals were zero (i.e. the dependent 
variable discretized algebraic equations were in balance) . 
The residual goes to zero as the number of sweeps goes to 
infinity. With numerical solutions, users must employ a 
variety of methods to determine problem convergence since 
residuals can not decrease to zero, as a practical matter. 

With the PHOENICS output (Appendix B) , this research 
utilized a local and a global analysis of the dependent 
variables to determine problem convergence. A local or 
"spot value" analysis can be accomplished with the PHOENICS 
output by designating cells to be monitored. Since the 
test facility field data referenced the exhaust gas 
temperature and velocity for a location in the augmentor 
tube/exhaust stack transition region, the computer model 
spot checked the dependent variable values at the same 
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location (the cell at IX=1, IY=5, and IX=22). From the 
initial program run to the final sweep (about 10,000), the 
spot values were closely monitored and compared. Local 
convergence at this cell occurred when little or no change 
in the dependent variable values was observed from sweep to 
sweep. Local or spot value convergence of dependent 
variables for one cell is a necessary but not sufficient 
condition for problem convergence. For this reason, the 
results at "ALL" cells were spot checked during the 
solution process. 

During the problem solution process, the PHOENICS code 
sums the dependent variable whole-field residuals. An 
examination of the whole-field sums provides a means of 
global investigation to determine problem convergence. If 
whole-field residuals decrease towards zero, with 
increasing iterations, problem convergence is suggested. 
However, if the residuals increase, with increasing 
iterations, problem divergence is occurring and some 
adjustment of the PHOENICS relaxation parameters is 
required. Global analysis is also accomplished by 
examination of the energy balance of the system. The 
dependent energy variable enthalpy (HI) residual is first 
divided by 1,000,000. This result is then divided by the 
absolute value of the sum of the enthalpy (HI) net sources. 
If the quotient decreases to less than about 0.05 (5 
percent error) , problem convergence is suggested. If the 
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quotient continues to increase, problem divergence is 
occuring and PHOENICS relaxation adjustment is again 
required . 

F. SUMMARY 

Although there was limited test facility field data to 
compare with numerical results, it is suggested that the 
computer modeling of the Denmark gas turbine test facility 
with PHOENICS code was successful. The field test data 
compared well with the numerical results. During the 
common sense examination, the conditions that existed in 
the five major regions of investigation made good sense and 
were indicative of what one might expect in a gas turbine 
test facility. The total mass flow continuity check 
indicated less than 0.003 percent error. Finally, local 
and global convergence was indicated. A system energy 
balance was within a 5 percent error. 
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VI. RECOMMENDATIONS 



A. INTRODUCTION 

It is believed the PHOENICS code is a useful tool for 
the solution of gas turbine flow simulation problems. 
However, since the Denmark test facility has not "FINELY" 
modeled the following computer model enhancements are 
suggested: BFC , PHOTON, alternate turbulence models, grid 

refinement, and baffle modeling. 

B. BFC 

Body fitted coordinates (BFC) are necessary to 
incorporate the cylindrical and curved surfaces of the jet 
engine, augmentor tube, and turning vanes. With the proper 
geometric shapes in the test facility model, mass flow in 
the augmentor tube inlet, the secondary air inlets, and the 
exhaust stack will be more representative of the actual 
test facility. 

C . PHOTON 

The implementation of the PHOENICS PHOTON option will 
greatly enhance numerical results investigation. A two or 
three dimensional representation of a dependent variable 
field is much more useful than a cell by cell examination 
of numbers. 
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D. ALTERNATE TURBULENCE MODELS 



One or more of the alternate PHOENICS turbulence models 
should be implemented and investigated. Since it was 
suggested that turbulence modeling is not an exact science, 
verifying the effect of the turbulence model on the 
solution is crucial. It may be that the turbulence model 
selected does not greatly effect the solution. 

E. GRID REFINEMENT 

Increasing the number of cells (grid refinement) in the 
system will enhance the problem solution process. During 
the slab-sweep iterative process, the dependent variable 
discretized algebraic equations will be solved in a "LESS 
ABRUPT" manner which will promote problem convergence. 

Grid refinement is also necessary to verify the solution is 
grid independent. 

F. BAFFLE MODELING 

Finally, by utilizing additional PHOENICS CONPOR 
commands, baffle material could be modeled into the 
primary/secondary air inlets and the exhaust stack. By 
experimenting with different values of porosity, the total 
mass flow rate of the computer model could more closely 
represent the total mass flow rate of the actual test 
facility. 
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APPENDIX A 



Q1 CODE 
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TALK= 


F >RUN( 1>1 ) 






* 






* 


* 






* 


* 






* 


* 


DENMARK 


JET-ENGINE TEST FACILITY COMPUTER MODEL 


* 


* 






* 


* 




NON-PHOENICS USER-DEFINED VARIABLES 


* 


* 






* 


* 


XLBL1 


THE FOLLOWING VARIABLES ARE USED TO DESCRIBE THE 


* 


* 


XLMID1 


PHYSICAL CHARACTERISTICS OF THE DENMARK TEST 


* 


* 


XLMID2 


FACILITY. THE FIRST LETTER OF THE VARIABLE 


* 


* 


XLMID3 


SPECIFIES THE DIRECTION IN THE COORDINATE 


* 


* 


XLBL2 


SYSTEM. 


* 


* 


XLBL3 


THE FOLLOWING CONVENTIONS ARE USED TO DESCRIBE 


* 


* 


XLEXT1 


A PARTICULAR ZONE OF THE MODEL: 


* 


* 


XLEXT2 


L LENGTH IN A COORDINATE DIRECTION 


* 


* 


XNCBL1 


MID MIDDLE OF SAY THE AUGMENTOR TUBE 


* 


* 


XNCMD1 


BL BOUNDARY LAYER 


* 


* 


XNCMD2 


EXT EXTERIOR ZONE OF STRUCTURE 


* 


* 


XNCMD3 


NC NUMBER OF CELLS IN A DIRECTION 


* 


* 


XNCBL2 


INT INTERNAL ZONE OF STRUCTURE 


* 


* 


XNCBL3 




* 


* 


XNCEX1 


EXAMPLE: XLBL1 THE LENGTH OF THE FIRST 


* 


* 


XNCEX2 


ZONE IN THE X DIRECTION. 


* 


* 


YLMID 




* 


* 


YLBL1 




* 


* 


YLBL2 


YNCEX2 THE NUMBER OF CELLS FOR THE 


* 


* 


YLEXT1 


SECOND EXTERNAL ZONE IN THE Y 


* 


* 


YLEXT2 


DIRECTION. 


* 


* 


YLBL3 




* 


* 


YLEXT3 


ZLINTS THE LENGTH OF THE FIFTH 


* 


* 


YLBL4 


INTERNAL ZONE IN THE Z 


* 


* 


YNCMID 


DIRECTION. 


* 


* 


YNCBL1 




* 


* 


YNCBL2 




* 


* 


YNCEX1 




* 


* 


YNCEX2 




* 


* 


YNCBL3 




* 


* 


YNCEX3 




* 


* 


YNCBL4 




* 


* 


ZLEXT1 




* 


* 


ZLINT1 




* 


* 


ZLINT2 




* 


* 


ZLINT3 




* 


* 


ZLINT<+ 




* 


* 


ZLINT5 




* 


* 


ZLINT6 




* 


* 


ZLEXT2 




* 


* 


ZNCEX1 




* 


* 


ZNCIN1 




* 


* 


ZNCIN2 




* 


* 


ZNCIN3 




* 


* 


ZNCIN4 




* 



* 


2NCIN5 




* 


* 


ZNCIN6 




* 


* 


ZNCEX2 




* 


* 






* 


* 






* 


* 


TMPJET 


TEMPERATURE OF JET EXHAUST (K) 


* 


* 


CPJET 


SPECIFIC HEAT OF JET EXHAUST <J/KG*K) 


* 


* 


PRSJET 


PRESSURE OF JET EXHAUST (N/M**2) 


* 


* 


GASCON 


GAS CONSTANT OF JET EXHAUST ( N*M/KG*K ) 


* 


* 


MASJET 


MASS FLOW RATE OF JET EXHAUST (KG/SEC) 


* 


* 


VELJET 


VELOCITY OF JET EXHAUST (M/SEC) 


* 


* 


RHOJET 


DENSITY OF JET EXHAUST (KG/M**3) 


* 


* 






* 


* 


TMPAMB 


AMBIENT AIR TEMPERATURE <K) 


* 


* 


PRSAMB 


AMBIENT AIR PRESSURE ( N/M**2 ) 


* 


* 


CPAMB 


SPECIFIC HEAT AMBIENT AIR ( J/KG*K ) 


* 


* 






* 


* 






* 


* 






* 


* 






* 


* 






* 


* 






* 


* 






* 



********************* *********************************************** 



* 

* 

* 

* 

* 

* 



VARIABLE DECLARATIONS 



PHOENICS ALLOWS 45 USER DEFINED REAL-NUMBER VARIABLES 



REAL ( XLBL1 ,XLMID1 ,XLMID2 ,XLMID3 ,XLBL2 >XLBL3 ,XLEXT1 ,XLEXT2 ) 

REAL( YLMID , YLBL1 , YLBL2 >YLEXT1 ,YLEXT2 , YLBL3 ,YLEXT3 ,YLBL4 ) 
REAL(ZLEXT1,ZLINT1,ZLINT2,ZLINT3>ZLINT4,ZLINT5,ZLINT6,ZLEXT2 ) 
REAUTMP JET, CPJET, PRSJET, GASCON, MASJET , VELJET , RHOJET ) 

RE AL ( TMPAMB , PRSAMB ,CPAN3 , DE LT , FACT ) 

REALtKEINJ,EPINJ,KEINA,EPINA ) 

* 

* PHOENICS ALLOWS 45 USER DEFINED INTEGER VARIABLES 

* 

INTEGER( XNCBL1 ,XNCMD1 ,XNCMD2 ,XNCMD3 ,XNCBL2 ,XNCBL3 ,XNCEX1 ,XNCEX2 ) 
INTEGER( YNCMID , YNCBL1 >YNCBL2 ,YNCEX1 ,YNCEX2 ,YNCBL3 ,YNCEX3 , YNCBL4 ) 
INTEGER(ZNCEX1,2NCIN1,ZNCIN2,ZNCIN3,2NCIN4,ZNCIN5,ZNCIN6,2NCEX2) 
* 

* 

* 



* 

* 

* 



* * 

* INITIALIZATION OF USER-DEFINED VARIABLES * 

* * 

V/ v V/ yy yy V/ \/ v V V V \/ vv s/ WV v vr w \/V\/ wyvvvvvvvv ^ V/ ^ V/ V \Z \/ V yv/ \y w \/ \/ \y w y w W V v \/ \/ \/ 

A A * A A A A^vAA AA A AA A^\^\ A A^\^A^\A AAA ^VA A A A A A A A aTA A A 

* * 
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* X DIRECTION LENGTHS (METERS) 


* 


# 


* 


XLBL1=0 .5 




XLMID1=0 .5 




XLMID2=0.5 




XLMID3=0 .5 




XLBL2=0 .5 




XLBL3-5 . 4 




XLEXT1=2.0 




* 


* 


* NUMBER OF CELLS IN X DIRECTION 


* 


* 


* 


XNCBL1=1 




XNCMD1=1 




XNCMD2=1 




XNCMD3=1 




XNCBL2=1 




XNCBL3=3 




XNCEX1=1 




* 


* 


* Y DIRECTION LENGTHS (METERS) 


X 


* 


X 


YLBL1-2 . 0 




YLMID=1.0 




YL3L2-2 . 0 




YLEXT1=3 . 85 




YLEXT2 =2 . IS 




YLBL3=6 . 0 




X 


X 


* NUMBER OF CELLS IN THE Y DIRECTION 


X 


X 


X 


YNCBL1=4 




YNCMID=2 




YNCBL2=4 




YNCEX1=2 




YNCEX2=1 




YNCBL3-3 




X 


X 


* Z DIRECTION LENGTHS (METERS) 


X 


X 


X 


ZLEXT1=2 . 0 




ZLINT1=23 . 5 




ZLINT2=3 . 2 




ZLINT3=1 . 5 




ZLINT4=6.0 




ZLINT5=16 . 3 




ZLINT6=3.4 




ZLEXT2=4.0 




X 


X 


* NUMBER OF CELLS IN THE Z DIRECTION 


X 


* 


X 


ZNCEX1=1 




ZNCIN1=5 




ZNCIN2=1 




ZNCIN3=2 
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ZNCIN4=6 

ZNCIN5=5 

ZNCIN6=6 

ZNCEX2=2 

* * 

* * 

* * 

U \rvi’ VU VVWVV\/\/V\/V\/VUWv/ \£ VV\/V/VV\/VV\/VV/N/VVUVU V/U V U\/U VVV/i/V/ V\/v/\/VVVVVVV U U V/ u 
A A*7>^A A 7WV AA^7W\ 7V7V AA A A7V A A A7V 7V 7V A AA/T A /P A A A AAA AAA A A A A A aTa A A A AA A A aTa AA AA A A A A A 

V/ W V V VZ V/ U W N* 1 w v/ V V \/ U U W U V V V/ V/ l/vx v/ v/ \/ W \/ V V W i/ VV/V^V/ y/ VV/ V/ V/ V V */ V V/ V V v/ V/ w W «/ W V 1/ y V V.U y/ v/ 
A AAA AAA A AA AA A A A A A A AAA AAA AAA A A A AAAAA A A AA AA A A AAA A AA AAA AAA A A A A7\A A A A AAA A 

* * 

* GROUP 1. RUN TITLE AND OTHER PRELIMINARIES * 

* * 

* * 

TEXTl DENMARK JET-ENGINE TEST FACILITY MODEL) 

* * 

* * 

* * 

WVVUWVVv/VV/v/ V v/v^v/VV^^VVVVW\/WV V/ VV/ v/ V/ V/ V yy vyyy \/V/V\/v/\/V/VVVV/\/S/^WVVVV/ V V/ V V/ 

"a^Paia A A^A A a A a a^a a A^P^A a a^P a ^Pa*a a ^>a*a a a A A a a/\a a a a a /va a a a'a ^aaT a a A^A A A A A A 

******************************************************************** 
* * 

* GROUP 2. TRANSCIENCE > TIME -STEP SPECIFICATION * 

* * 

* DEFAULT TAKEN: STEADY STATE * 

* * 

* * 

* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 3. X-DIRECTION GRID SPECIFICATION * 

* * 

******************************************************************** 
* * 

* NX = THE NUMBER OF CELLS IN THE X DIRECTION (INTEGER) * 

* XULAST = DISTANCE IN THE X DIRECTION (METERS) * 

* * 
NX = XNCBL1 +XNCMD 1 +XNCMD 2 +XNCMD3 +XNCBL 2+XNCBL3 +XNCEX1 

* 

XULAST =XL3L1 +XLMID1+XLMID2+XLMID3+XLBL2 + XLBL3+XLEXT1 

* * 

* THE METHOD-OF-PAIRS HAS BEEN USED TO DEFINE THE NUMBER * 

* AND SIZE OF CELLS. THE XFRAC >YFRAC > AND ZFRAC SPECIFY THE X,Y>* 

* AND Z COORDINATES RESPECTIVELY. THE ODD INDEX SPECIFIES THE * 

* NUMBER OF CELLS IN A DIRECTION AND THE EVEN INDEX SPECIFIES THE* 

* SIZE OF EACH CELL. THE NEGATIVE QUANTITY FOR THE FIRST ODD * 

* INDEX IN EACH GROUP IS THE RESERVED WORD WHICH SPECIFIES * 

* "METHOD-OF-PAIRS*' * 

* * 

XFRAC ( 1 ) = -XNCBLl \ XFRAC( 2 )=XLBLl/( XNCBL1*XULAST ) 

XFRAC ( 3 ) =XNCMD1 > XFRAC ( <t ) =XLMI D 1/ ( XNCMD I*XULA$T ) 

XFRAC ( 5 )=XNCMD2 * XFRAC ( 6 ) =XLMID2/( XNCMD2*XULAST ) 

XFRAC ( 7 )=XNCMD3 > XFRAC ( 8 )=XLMID3/( XNCMD3*XULAST ) 

XFRAC ( 9 ) =XNCBL2 > XFRAC (10 ) =XLBL2/( XNCBL2*XULAST ) 

XFRAC ( 11 )=XNCBL3 > XFRACt 12 )=XLBL3/( XNCBL3*XULAST ) 
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XFRAC ( 13 ) =XNCEX1 > XFRAC ( 14 ) =XLEXTl/( XNCEX1*XULAST ) 

* * 



* * 

* * 

■VV V V V V V y/\y \/U\/V \/V\/ VV V V W V V V \/ V V V vr V vv^v V V V^V V V vwv/w^Vwv/V\yv vUwv/Vwv/ w 

A ^ A A A A A A A A AA^A A A A A^A A A^AA ^AA A^^^xA A^AA A ^AA A A^A A^AA A A A^A A A A a7\ A^^AA A7\ A A A A^\ a7\ A a7\ A 



* * 



* GROUP 4. Y-DIRECTION GRID SPECIFICATION * 




* 

* NY = THE NUMBER OF CELLS IN THE Y DIRECTION (INTEGER) 

* YVLAST = DISTANCE IN THE Y DIRECTION (METERS) 

* 

NY=YNCMID+YNCBL1+YNCBL2+YNCEX1+YNCEX2+YNCBL3 

YVLAST=YLMID+YLBL1+YLBL2+YLEXT1+YLEXT2+YLBL3 

* 

* METHOD OF PAIRS 

* 



YFRACt 1 ) = -YNCBLl $ 
YFRAC ( 3 )=YNCMID > 
YFRAC( 5 ) = YNCBL2> 
YFRAC( 7 )=YNCEX1> 
YFRAC ( 9 )=YNCEX2 > 
YFRAC( 11 )=YNCBL3) 
* 



YFRAC ( 2 )=YLBLl/( YNCBL1*YVLAST ) 
YFRAC ( 4 )=YLMID/( YNCMID*YVLAST ) 
YFRAC ( 6 )=YLBL2/( YNCBL2*YVLAST ) 
YFRAC ( 8 )=YLEXTl/( YNCEX1*YVLAST ) 
YFRAC( 10 )=YLEXT2/( YNCEX2*YVLAST ) 
YFRAC ( 12 )=YLBL3/( YNCBL3*YVLAST ) 



* 

* 

* 

* 



* 

* 

* 



* 



* * 
******************************************************************** 
******************************************************************** 



GROUP S. 2-DIRECTION GRID SPECIFICATION 



******************************************************************** 
* * 

* N2 = THE NUMBER OF CELLS IN THE 2 DIRECTION (INTEGER) * 

* 2WLAST = DISTANCE IN THE 2 DIRECTION (METERS) * 

* * 
NZ=ZNCEX1+ZNCIN1+ZNCIN2+ZNCIN3+ZNCIN4+ZNCIN5+ZNCIN6 +2NCEX2 

2WLAST =2LEXT1+2LINT1 +2LINT2 + 2LINT3+2LINT4 + ZLINT5+2LINT6 +2LEXT2 



2FRAC( 1 ) = -ZNCEXl > 
2FRAC ( 3 )=ZNCIN1 > 
2FRAC ( 5 )=ZNCIN2 > 
ZFRAC( 7)=ZNCIN3> 
2FRAC ( 9 ) =ZNCIN4 > 
2FRAC ( 11 )=2NCIN5 > 
ZFRAC( 13 )=ZNCIN6; 
2FRAC ( 15 )=ZNCEX2 > 
* 



ZFRAC( 2 )=ZLEXTl/( ZNCEX1*ZWLAST ) 
2FRAC ( 4 )=ZLINTl/( ZNCIN1*ZWLAST ) 
ZFRAC ( 6 )=ZLINT2/( ZNCIN2*ZWLAST ) 
ZFRAC( 8 )=ZLINT3/( ZNCIN3*ZHLAST ) 
ZFRAC ( 10 )=ZLINT4/(ZNCIN4*ZWLAST) 
ZFRAC ( 12 ) =2LINT5/( ZNCIN5*ZWLAST ) 
ZFRAC ( 14) =ZLINT6/(2NCIN6*2WLAST ) 
2FRAC( 16 )=2LEXT2/( ZNCEX2*ZWLAST ) 



* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 6. BODY-FITTED COORDINATES OR GRID DIRECTION * 
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* ( NOT IN USE ) * 

Vtf^/vv^x/N/N/uyyyyyyy vvv w vv vw vv^vvvvvvvv u wi/vvvw^wwx/N/w/x/v v vvvvvvMv viwy 
Tv A M ^T^V^V^V A^T A AAA A A A A A A A AAA A A A A A^AA^ A A Aa A A^VA a!aa A^V A*AaTa A A A A A A A^T A 

* * 

* * 

* * 

* * 

* GROUP 7. VARIABLES STORED, SOLVED AND NAMED * 

* * 

* * 

SOLVE ( P1>H1,U1,V1>W1 ) 

SOLUTNt P1,Y,Y,Y,N>N,N) 

STORE ( RH01 ,TMP1 ) 

* * 

* * 

* * 

******************************************************************** 
* * 

* GROUP 8. TERMS (IN DIFFERENTIAL EQUATION) AND DEVICES * 

* * 

WVUvW VWVVVVVVVVVVMVtfWtfVVVV^WWV.WWVt/V^WV^ VfV WWVVV N/Wk/V t/VtfMV VWV/N/ 
A A^V A A A A A A A A A A A A A AAA A A A A ^»^VA^V A A^V^V^VaTA A^V A A A^V A A A AA^VAAA A A A^V A^V A A A^V A 

* * 

* DEFAULT TAKEN * 

* * 

* * 

* * 

******* ************************************************************* 
******************************************************************** 
* * 

* GROUP 9. PROPERTIES OF THE MEDIUM (OR MEDIA) * 

* * 

********************************* x ************** ******************** 
* * 

TMP1A=TINY 
TMP AMB = 295. 

PRSAMB=0 . 

CPAMB=1004. 

TMP IB =1 . 0/CP AMB 
TMP1=GRND2 
TMP JET =2000 . 

CPJET=1004. 

PRSJET=0 . 

GASC0N=287. 

RH01B=1. /GASCON 
MASJET=100 . 

PRESSQ=I01325. 

RH0JET = ( PRESSO + PRSJET )/( GASCON*TMP JET ) 

VELJET=MA3JET/RH0JET 
RH01=GRND5 
ENUL-1. 0E-05 
DRH1DP=GRN05 
TURMCD( KEMODL ) 
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STORE (VIST) 

KEIMJ=0.5*( 0.005*VELJET)**2 
E PINJ = 0 . 09*KEINJ**1 .5/(0. 05*0 . 5 ) 
KEINA=0 . 5*( 0 . 005*1 . 0 )**2 
EPINA r 0 . 09*KEINA**1 . 5/( 0.1) 




* GROUP 10. INTERPHASE -TRANSFER PROCESSES AND PROPERTIES * 



* * 



VW*/\/\/VtfVV VV ’^W V VV\/ VV l VV\/\/\/VWVV \/VVVVVU\/VVVVVV\/\/V\/\/tfVV\/VVV\/V\f\fWV WVV\/ 
* * * a! a a a aaT a a a a a a a a a^Pa^a a a a a a a ^Pa a a a aAa? a"a^Paa^1aa ^Paa^ a^P^P ^Pa 

* DEFAULT TAKEN * 

* * 



V VVVvVvVVVSf.VVVVVVy VV’W V vv vVV vvVv^Vvv/wvWVV v vvvvvWvuVVVVVvvW vVV/VV 

A A^AA^ A A A a A A AiAA^ A A A^P A A A ^Pa AAA A A AAA AAA A A^P A^AA^ A AiAiA A A A A A A A ^Ta A^P/Pa A A A A A A 7* ^P»P »\ 

\/VV\/i*v#V/VWW\/V>/V'WVVW *j' WVW.V v V/ \/\/V \^WVU%/VWUWVVW V VWWW \/ V V V W V V V v vV 

AATAA AA AA A AA A /P/PaA 7\ AA /P/Pa A AA A A7NA AA A A AAA AATPa A AAVPa 7P A 7P7P A AAAA AAA AA AA A A /Pa/P A 

* * 

* GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS * 

* * 

VW^yN/WVWWVVV/V WWVV^W V W W V W ^ V X/WWVVWWVV^WWWVV/VV/V X/VUW W VV/VV 

A A A A A A A A AAA A A a AAA A AAA A a A A A A A A A AAA AAA A A AA A A A AAA AA A A A A A AAAA A AAAA A A AAA A 

* * 

FIINIT( HI )=CPAMB*TMPAMB 
FIINIT(TMP1)=TMPAMB 
FIINIT(DEN1)=1.0 
FIINIT( PI ) = 0 . 0 
FIINIT(U1 )=0.0 
FIINIT( VI )=0 . 0 
FIINIT( HI )=0 . 0 
FIINIT( KE )=KEINA 
FIINIT(EP)=EPINA 

* * 



* 



* 



* * 

* BLOCKING OUT THE MASS OF THE JET ENGINE * 

* * 
CONPOR(0.0,CELL,1,1>5,6,6,7) 

* * 



* NALL, FLOOR, AND CEILING BOUNDARY CONDITIONS. RESTRICTING FLOW* 

* ACROSS A BOUNDARY BY USING CONPOR VALUES FROM 0.0 (NO FLOW) TO* 

* 1.0 (MAX FLOW). BOUNDARY FRICTION TO FLOW IS ALSO IMPLEMENTED* 

* USING A NEGATIVE BOUNDARY COORDINATE. THIS CONVENTION CAUSES * 

* FLOW PARALLEL TO THAT BOUNDARY FACE TO HAVE FRICTION. BOUNDARY* 



* CONDITIONS HILL BE DESCRIBED BY COMMENT AND PLANE. * 

* * 

* * 

* TEST FACILITY FOUNDATION AND GROUND OUTSIDE STRUCTURE * 

* * 



CCNP0R( 0 . 0 , SOUTH , 1 ,NX ,1 , 1 , 1 ,NZ ) 

PATCH( FL0R, SHALL, 1, NX, 1,1,1,NZ, 1,1) 

COVAL ( FL0R,U1 >GRND2 ,0.0) 

C0VAL( FLOR ,N1 ,GRND2 ,0.0) 

COVAL ( FLOR ,KE ,GRND2 ,GRND2 ) 

C0VAL( FL0R,EP ,GRND2,GRND2 ) 

* * 
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JET-ENGINE TEST BAY 



* 

* 

* 

* 



BAY SIDE WALL !Y-Z) 



CONPOR! 0.0, WEST, 9, 9, 1,12, 2, 9) 
* 



* BAY CEILING !X-Z) 

* 



CONPOR( 0 . 0 , SOUTH , 1 ,8 , 13 , 13 , 2 , 9 ) 

* 

* BAY WALL REAR (X-Y) 2 SECTIONS 

* 

CONPOR( 0.0, LOW, 6, 8, 1,12,-10,10) 

CONPOR( 0.0, LOW ,1,5,11,12,-10,10) 

* 

* EXHAUST DUCT INLET TUBE 

* 

CONPOR( 0 . 0 , NORTH ,1,2,3,3,10,14) 

PATCH ( TUB1 ,NWALL ,1,2,3,5,10,14,1,1) 
COVAL( TUB1 >U1 ,GRND2 ,0.0) 
COVAL(TUB1,W1,GRND2,0.0) 

COVAU TUB1 ,KE ,GRND2 ,GRND2 ) 

COVAL! TUB1 ,EP ,GRND2 ,GRND2 ) 

CCMPOR( 0 . 0 , SOUTH ,1,2,8,8,10,14) 

PATCH! TU32 ,SWALL , 1,2,8,3,10,14,1,1) 
COVAL ( TUB 2 ,U1 ,GRND2 , 0.0) 

COVAL( TUB 2 ,W1 ,GRND2 ,0 . 0 ) 

COVAL! TUB 2 ,KE ,GRND2,GRND2 ) 

COVAL! TUB2 , E P , GRND 2 , GRND 2 ) 

CONPOR! 0.0, WEST, 3, 3, 4, 7, 10, 14 ) 

PATCH! TUBS, WWALL, 3, 3, 4, 7, 10, 14, 1,1 ) 
COVAL ! TUBS , VI , GRND 2,0.0) 

COVAL! TU33 ,W1 ,GRND2 ,0 . 0 ) 

COVAL! TUB3 ,KE ,GRND2 , GRND 2 ) 

COVAL ( TUB3 , EP , GRND 2 , GRND 2 ) 

CONPOR! 0 . 0 , LOW , 1 ,5 , 1 , 3 , 10 , 10 ) 

CONPOR! 0.0, LOW, 1,5, 8, 10, 10,10 ) 

CONPOR! 0.0, LOW, 3 ,5 ,4 ,7,10,10) 

* 

* SECONDARY-AIR INLET WALL tX-Y) 

* 



CONPOR! 0 .0 ,LCW,6 ,7,1,10,12,12 ) 

PATCH! SAIW , LWALL ,6,7,1,10,12,12,1,1 ) 
COVAL! SAIW, U1,GRND2, 0.0) 

COVAL (SAIW, VI, GRND2, 0. 0 ) 

COVAL! SAIW, KE,GRND2,GRND2 ) 

COVAL! SAIW , EP , GRND 2 , GRND 2 ) 

* 



* 

* 

* 

* 

* 

* 

* 

* 

* 

* 



* 

* 

* 



* 

* 

* 



* SECONDARY-AIR INLET AND PARTIAL EXHAUST DUCT CEILING ( X-Z ) * 



* * 
CONPOR! 0.0, SOUTH, 1,7, 11, 11, 10, 12 ) 

PATCH! SAIWPE ,SWALL ,1,7,11,11,10,12,1,1 ) 

COVAL ( SAIWPE ,U1 ,GRND2 ,0 . 0 ) 

COVAL ! SAIWPE ,W1 ,GRND2 , 0.0) 

COVAL! SAIWPE ,KE ,GRND2 ,GRND2 ) 
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COVAL! SAIWPE,EP,GRND2,GRND2) 

* 

* EXHAUST DUCT CEILING ( X-Z ) 

* 

CONPOR! 0.0, SOUTH, 1,5, 11, 11, 13, 20 ) 

PATCH! EDC, SWALL, 1,5, 11, 11, 13, 20, 1,1) 

COVAL! EDC ,U1 ,GRND2 ,0 . 0 ) 

COVAL! EDC ,W1 ,GRND2 ,0.0) 

COVAL! EDC ,KE ,GRND2 ,GRND2 ) 

COVAL( EDC,EP,GRND2,GRND2) 

* 

* EXHAUST-DUCT WALL (Y-Z) 

* 

CONPOR! 0.0, WE ST, 6, 6, 1,10, 13, 20) 

PATCH! EDWXZ , WWALL ,6 ,6,1,10,13,20,1,1) 

COVAL! EDWXZ , VI ,GRND2 ,0.0) 

COVAL( EDWXZ ,W1 ,GRND2 ,0.0) 

COVAU EDWXZ, KE ,GRND2 ,GRND2 ) 

C OVA L ( E DWXZ , E P , GRND 2 , GRN D 2 ) 

* 

* EXHAUST-STACK WALLS (X-Y) TWO SECTIONS 

* 

CONPOR! 0.0, HIGH, 1,5, 11, 13, 20, 20) 

PATCH! ESWXY1 ,HWALL ,1,5,11,13,20,20,1,1) 

COVAL! ESWXY1 ,U1 ,GRND2 , 0.0) 

COVAL ( ESWXY 1 , VI , GRND2 ,0.0) 

COVAL! ESWXY1,KE , GRND 2 ,GRND2 ) 

COVA L ( E SWXY 1 , E P , GRND 2 ,.GRND 2 ) 

CONPOR! 0.0, LOW, 1,5, 1,13, 27, 27) 

PATCH! ESWXY2,LWALL, 1,5, 1,13, 27, 27, 1,1) 

COVAL! ESWXY 2 ,U1 ,GRND2 ,0.0) 

COVAL! ESWXY 2 ,V1 ,GRND2 ,0 . 0 ) 

COVAL! ESWXY2,KE ,GRND2 ,GRND2 ) 

COVAL! ESWXY2 ,EP , GRND 2, GRND 2 ) 

* 

* EXHAUST-STACK WALL (Y-Z) 

* 

CONPOR! 0.0, WEST, 6, 6, 1,13, 21, 26 ) 

PATCH! ESWXZ, WWALL, 6, 6, 1,13, 21, 26, 1,1) 

COVAL ! ESWXZ , VI ,GRND2 ,0.0) 

COVAL ! ESWXZ ,Wi ,GRND2 , 0.0) 

COVAL! ESWXZ, KE ,GRND2 ,GRND2 ) 

COVA L ! E SWXZ , E P , GRND 2 , GRND 2 ) 

* 

* TURNING VANES 



* 

* 

* 



* 

* 

* 



* 

* 

* 



* 

* 

* 



* 

* 



* * 
CONPOR! 0.0, SOUTH, 1,5, -7, 7, 21, 22) 

CONPOR (0.0, SOUTH ,1,5,-3,3,23,24) 

CONPOR! 0.0, LOW, 1,5, 7, 10, -23, 23 ) 

CONPOR ( 0 . 0 , LOW , 1 ,5 , 3 , 6 , - 25 , 25 ) 

**********^***x****************^*****~**-**********~* **************** 

***************** *-*X-***^*-X*^BH***"*********^*XX XXX XXX****XX***X***"***"X 



* * 

* GROUP 12. CONVECTION AND DIFFUSION ADJUSTMENTS * 

* * 
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DEFAULT TAKEN 




* 

* 



* GROUP 13. BOUNDARY CONDITIONS AND SPECIAL SOURCES 




JET-ENGINE MASS SOURCE 



* 

* 



PATCH( JETIN, HIGH, 1,1, 5, 6, 5, 5, 1,1 ) 

COVAL ( JETIN , P 1 , FIXFLU , -MASJET*1 . ) 

COVAL( JETIN , HI ,ONLYMS ,SAME ) 

PATCH! JETOUT, LOW, 1,1, 5, 6, 8, 8, 1,1) 

COVAL ( JETOUT , PI , FIXFLU, MASJET*1 . J 
COVAL ( JETOUT ,W1 ,ONLYMS ,VE LJET ) 

COVAL ( JETOUT ,H1 ,QNLYMS,CPJET*TMP JET ) 

COVAL! JETOUT ,KE ,ONLYMS ,KEINJ ) 

CO VA L ( J E TOUT , E P , ON L YMS , E P I N J ) 

* * 

* ATMOSPHERIC BOUNDARY CONDITIONS * 

* * 

* LAST CELL SURFACE IN X DIRECTION * 

* * 
PATCH! XSKY >EAST , NX, NX, 1, NY ,1,NZ,1,1) 

COVAL ( XSKY , PI, 0. 1,0.0 ) 

COVA L ( XSKY , H 1 , ONL YMS , C PAMB*TMP AMB ) 

C 0 V A L ( XSK Y , K E , ON L YMS , K E I NA J 
COVA L ( XSKY , E P , ONL YMS , E P IN A J 

* * 

* LAST CELL SURFACE IN THE Y DIRECTION * 

* * 
PATCH( YSKY , NORTH , 1 ,NX ,NY ,NY , 1 ,N2, 1 , 1 ) 

COVAL ( YSKY , PI, 0.1, 0.0) 

COVAL ( YSKY , HI ,ONLYMS , CPAMB*TMP AMB ) 

COVAL! YSKY ,KE ,ONLYMS ,KEINA ) 

COVAL! YSKY, EP ,ONLYMS ,EPINA ) 

* * 



PATCH! ZSKY1 , LOW , 1 ,NX ,1, NY, 1,1, 1,1) 

COVAL! 2SKY1, PI, 0. 1,0.0 J 

COVAL! 2SKY1 ,H1 ,ONLYMS ,CPAMB*TMPAMB ) 

COVAL! ZSKY1 ,KE ,ONLYMS,KEINA ) 

COVAL ( ZSKY1 , E P , ONL YMS , E PINA ) 

* * 

* LAST CELL SURFACE IN THE 2 DIRECTION * 

* * 
PATCH! ZSKY2, HIGH, 1, NX, 1, NY, NZ,NZ, 1,1) 

COVAL! ZSKY2, PI, 0. 1, 0. 0 ) 

COVA L ( 2SKY 2 , H 1 , ONL YMS , CPAM3*TMPAMB ) 

COVAL ! 2SKY2 ,KE ,ONLYMS ,KEINA ) 



* FIRST CELL SURFACE IN THE Z DIRECTION 



* 



* 
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COVAL( ZSKY2 ,EP ,ONLYMS,EPINA ) 




* DEFAULT TAKEN * 



* * 

V V \/ A/ A/ V VV V! V v v v UV \/ W V/ V ^V/VVWVVWWVV V VWWVvyvv/vyVWV WWWW v v v V \y \/ \/ \# V V 
« A A /T7\« AAA A A ^ A * *7\A7\ 7\ 7\ /V 7wT* /V 7 wT* AA^. A A! ■JT A A /V R R « AATAAA A AA AA AA aTaA A A A A 

W-VV vWx/^Vwv WWW WVV’WVVVAZVlfV VV V VW^ VLV_ W VVV UV\/WVV\/Vv V V W V VW-W-WWWVTWNZ 
A A A A A A A AA A A A A A A A A AAA AA A A AAAAA A AAA A A A A A AAA A A A A A A 7V A A A A A A AAA A A A AAA A AAAA 



* 

* GROUP 15. TERMINATION OF SWEEPS 

* 



* 

* 

* 



* * 

RESTRT(ALL) 

FSWEEP=1 

LSWEEP=2 

* * 



* 



* 



* * 

^■WVWW W V v v vv w vw v viviw: ^^ VV ^V W^ W W V V V VV VV VV VV vV VVVV VV VVVUVW V WV 

A R A K IUI A W A R A A ^ A K A A 7 TX 7 TKR H Rif A H K JinAAJlJl K’W AAARjl 

************************************************************* ** ***** 



* 

* GROUP 16. TERMINATION OF ITERATIONS 

* 



* 

* 

* 



V v w v\/\yvvv\/v vw \/ \z vv ^w vl w vv v v v uw v v^vvvvw w v v:\zvvvw v ^ v v vv v w v \/ v v v V w 

A A A A A A A A A A A A AAA A A A AAA AAA A A A A A A A A A A A A A A A A A A A A AA AAA A A AAA AAA A A A A A A A A AAA 



* 



* 



* DEFAULT TAKEN 

* 



* 

* 



******************************************************************** 



******************************************************************** 
* * 

* GROUP 17. UNDERRELAXATION DEVICES * 

* * 

\/ VV V V V V VVVWV^/V^ V W.^VV V WWW\/WVWV\/VV \/ W V WWVWVVW \/V V\/WVVV\/W VV 
A AAA A A A AAA A A A A^V^\ A aAAT A A A A A A AAA A/\ A A 7\A^\AAA A A A A A A A A A A^V A A A A A A AA AA A A A A A 



* * 
DELT =£0 . /( NZ*0 . 2*VELJET ) 

FACT =0 . 25 

SHIFTED TO FACT=0.5 AT 1550 SWEEPS 
FACT = 0 .5 

RELAX(P1,LINRLX,0.3) 

RE LAX( U1 , FALSDT , FACT*DELT ) 

RE LAX( VI , FALSDT , FACT*DELT J 
RE LAX( W1 , FALSDT > FACT*DELT ) 

RE LAX( KE , FALSDT , FACT*DE LT ) 

RE LAX ( EP, FALSDT, FACT*DELT ) 

REMOVED THE 2.5 FROM HI RELAX AT 5000 SWEEPS 
RE LAX ( HI , FALSDT, FACT*DELT*1 . 0 ) 

* * 



* * 

* * 

WW^ w Wx/A^N/V^v v Vy/v/x/vwWK/vv/VK/ v/vr wVN/VWVSfV WVl/\/\yVVVV\/V/V\/W\/V\/\/V\/V.VWW 
T' A A A A A A A A A A A AAA A A^AAAA A A AAA A AAA A AAA A A AAA AAA A A A AAA A ATT /\ A A A A A A A A A A A A A A A 
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U V V W V V V/ \/ V w v v yyvvvwvv W V WWW VVVVVX/ VVW\/\/V WV\/WS/VVVVVVVW\/ V V/ S/ W V 

A A. A A A A 7V 7VA 7\ 7\ 7\ A 7T7\A A A A AAA A A A *" A AA A A AA A A A A A A A A A A A A AAA AAA AA AAA A A A A A A AAA A 

* * 

* GROUP 18. LIMITS ON VARIABLES OR INCREMENTS TO THEM * 

* * 

WW¥V¥VV¥^ VVVV\/\/ WWVWWV\/\/ V V V V ^VVIVVVU V WVVVW^/VVVVWVWN/VVWWWt/S/VUV w 
A A /\ A A A A A A A7\ A A A A A A AA A A 7w\ A A A A A A A A A A A AA7\ A A A7\ A AAAA A A AAAA 7\A AAAA A AA AA 7\ A A7\ A 

* * 

VARMIN(P1)=-2*PRESS0 

VARMIN( HI )=CPAMB*TMPAKB 

VARMIN(U1)=-VELJET 

VARMINf VI ) = -VE LJET 

VARMIN(W1)=-VELJET 

VARMIN( TMP1 )=TMPAMB 

VARMIN( RH01 3=0 . 1 

VARMAX( P1)=2*PRESS0 

VARMAX( HI )=CPJET*TMPJET 

VARMAXt U1 )=VELJET 

VARMAX( VI )=VELJET 

VARMAX( HI )=VELJET 

VARMAX( TMP1 J=TMPJET 

VARMAX( RH01 )=2 . 0 

* * 

* * 

* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 19. DATA COMMUNICATED BY SATELLITE TO GROUND * 

* * 

******************************************************************** 
* * 

* DEFAULT TAKEN * 

* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 20. PRELIMINARY PRINTOUT * 

* * 

******************************************************************** 
* * 

* * 

* DEFAULT TAKEN * 

* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 21. PRINTOUT OF VARIABLES * 

* * 

******************************************************************** 
* * 

* DEFAULT TAKEN * 

* * 

******************************************************************** 
******************************************************************** 
* * 

* GROUP 22. SPOT-VALUE PRINTOUT * 
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IXM0N=1 

IYM0N=5 

IZMON=22 



TSTSHP=10 



* 

* 

* 



* 

* 

* 



******** XXXX ******************************************* * ************ 



PATCH! PL0T5, PROFIL, 1,1, 5, 5, 1,NZ, 1,1) 

PLOT( PLOTS, PI, -3000. 0,2500.0 J 
PATCH ( PL0T3,PR0FIL, 1,1, 5,5,1, NZ, 1,1 ) 

PL0T( PL0T3 ,U1 , -20 . 0 ,50 . 0 ) 

PATCH! PL0T4,PR0FIL,1,1,5,5,1,NZ,1,1) 

PLOT! PL0T4 , VI , -40 . 0 , 60 . 0 ) 

PATCH! PLOT2 ,PROFIL,l ,1,5,5,1,NZ,1,1) 

PLOT! PLOT2 ,W1 ,-15.0 ,500 . 0 ) 

PATCH! PL0T1,PR0FIL, 1,1, 5,5,1 ,NZ, 1,1) 

PLOT! PLOT1 , HI, 2750 00. 0,2100000.0) 

PATCH! PLOT6, PROFIL, 1,1, 1,10, 20, 20, 1,1) 

PLOT! PL0T6, HI, 300000.0, 1000000. 0 ) 

PATCH! PLOT7 ,PROFIL ,1,1,1,10,20,20,1,1 ) 

PLOT! PLOT7 ,W1 , 20 . 0 , 1 75 . 0 ) 

PATCH! PL0T8 , PROFIL, 1,5,5,5,20,20,1,1 ) 

PLOT! PL0T8, HI, 300000. 0,1000000. 0 ) 

PATCH! PL0T9, PROFIL, 1,5, 5, 5, 20, 20, 1,1) 

PLOT! P LOT 9 ,W1 , 20 . 0 , 1 75 . 0 ) 

PATCH! PLOTIO, PROFIL, 1,1, 1,10, 16, 16, 1, 1 ) 

PLOT! PLCTIO, HI, 300000. 0,1000000.0) 

PATCH! P LOTH, PROFIL, 1,1, 1,10, 16, 16, 1,1) 

PLOT! PLOT1 1 ,W1 , 20 . 0 ,175 . 0 ) 

PATCH! PLOT12, PROFIL, 1,5, 5, 5, 16, 16, 1,1) 

PLOT ( PLOT 12 ,H1 ,300000.0,1000000.0) 

PATCH! PLOT13, PROFIL, 1,5, 5, 5, 16, 16, 1,1) 

PLOT! PL0T13,W1, 20. 0,175.0) 

PATCH! PL0T14, PROFIL, 1,1, 13, 13, 21, 26, 1,1) 

PLOT! PL0T14 ,H1 ,500000 . 0 , 600000 . 0 ) 

PATCH! PLOT15, PROFIL, 1,1, 13, 13, 21, 26, 1,1) 

PLOT! PL0T15, VI, 10.0,45.0 ) 

PATCH! PLOT16, PROFIL, 1,5, 13, 13, 23, 23, 1,1) 

PLOT! PLOT16, HI, 500000. 0,600000. 0 ) 

PATCH! PLOT17, PROFIL, 1,5, 13, 13, 23, 23, 1,1) 

PLOT! PLOT 17, VI, 10. 0,45.0) 

NUMCLS=10*NXPRIN 
NPLT = 10 
ITABL=3 
IPROF=3 

LUPR3=6 

* * 

* * 

******************************************************************** 
************************************** * ***************** * * x x x x ****** 
* * 

* GROUP 24. DUMPS FOR RESTARTS * 

* * 




* 

* 



* GROUP 23. FIELD PRINTOUT AND PLOT CONTROL 




* 



* 




STOP 
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APPENDIX B 



EARTH OUTPUT 
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FIELD VALUES OF KE 

IY= 16 3.816E-01. 3.80SE-01 3.782E-01 3.749E-01 3.704E-01 3.579E-01 3.257E-01 2.829E-01 2.289E-01 
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FIELD VALUES OF KE 
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Figure C.l • Denmark Jet-Engine Test Facility 
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Figure C.2. Staggered-Grid Cell Configuration 
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Figure C.3. Cartesian Cell Configuration 
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HIGH « LOW 




GROUP TITLE 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 



RUN TITLE AND OTHER PRELIMINARIES 
TRANSIENCE: TIME-STEP SPECI FICATI ONS 
X-DIRECTION GRID SPECI FICATI ON 
Y-DIRECTION GRID SPECIFICATIONS 
Z-DIRECTION GRID SPECI FICATI ONS 
BODY-FITTED COORDINATES OR GRID DISTORTION 
VARIABLES STORED, SOLVED AND NAMED 
TERMS (IN DIFFERENTIAL EQUATIONS) AND DEVICES 
PROPERTIES OF THE MEDIUM (OR MEDIA) 
INTERPHASE-TRANSFER PROCESSES AND PROPERTIES 
INITIALIZATION OF VARIABLE OR POROSITY FIELDS 
CONVECTION AND DIFFUSION ADJUSTMENTS 
BOUNDARY CONDITIONS AND SPECIAL SOURCES 
DOWNSTREAM PRESSURE FOR PARAB=T 
TERMINATION OF SWEEPS 
TERMINATION OF ITERATIONS 
UNDERRELAXATION DEVICES 

LIMITS ON VARIABLES OR INCREMENTS TO THEM 

DATA COMMUNICATED BY SATELLITE TO GROUND 

PRELIMINARY PRINTOUT 

PRINTOUT OF VARIABLES 

SPOT-VALUE PRINTOUT 

FIELD PRINTOUT AND PLOT CONTROL 

DUMPS FOR RESTARTS 



Figure C.4. PHOENICS Data Group Catagories 
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Figure C.5. CONPOR Porosity And Blockage 
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FILE: C0PYQ1 EXEC 



A 1 



STRACE 

-START 

STYPE WHAT ARE THE FILENAME AND FILETYPE OF THE Q1 FILE YOU WANT 
STYPE TO STORE ON MVS? USE ENTER TO EXIT. 

SREAD VAR SQ1 SQ2 

SIF /SQ1 EQ / SEXIT 
STATE SQ1 SQ2 A 
& I F SRC EQ 0 SGOTO -NEXT 

STYPE FILE SQ1 SQ2 A DOES NOT EXIST. USE ENTER TO EXIT. 

SREAD VAR S 
SIF /S EQ / SEXIT 
SGOTO -START 
-NEXT 

STATE COPYQ1 JCL A 

SIF SRC EQ 0 SGOTO -COPY 

STYPE TEMPLATE FILE COPYQ1 JCL A DOES NOT EXIST ; 

STYPE SEE LAB MANAGER. 

SEXIT -2 

-COPY 

SERROR SEXIT -1 

COPYFILE COPYQ1 JCL A TEMPI JCL A ( REPL 
S = SCONCAT OF /YOURQ1 FI L ENAME/ SQ1 / 

SSTACK TOP 

SSTACK CL OCATE/YOURQ1 FI L ENAME/ 

SSTACK CHANGE S 1 1 
SSTACK CLOCATE/SYSUT1/ 

SSTACK GET SQ1 SQ2 A 
SSTACK TOP 
S STACK FILE 
XEDIT TEMPI JCL A 
EXEC SUBMIT TEMPI JCL A 
CONWAIT 
DESBUF 
SEXIT 



FILE: COPYQ1 JCL A1 



//CXHUSH JOB (4555,9999), ' COPYTMVS » , CL ASS = A 
// EXEC PGM= I EBGENER 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD DUMMY 

//SYSUT2 DD DSN=MSS . S4555 . PHOENICS . SRC (YOURQ1 FILENAME ) , DISP=SHR 
//SYSUT1 DD x 
/x 
// 



Figure C.6 COPYQ1 EXEC And JCL Files 



93 



FILE: RUNSAT EXEC A1 



&TRACE 

-START 

&TYPE WHAT IS THE FILENAME OF THE Q1 FILE YOU WANT TO RUN ( SATEL ITE) ? 
&TYPE USE ENTER TO EXIT. 

&READ VAR &Q1 

& I F /&Q1 EQ / &EXIT 
X STATE SQ1 DATA A 
x SIF SRC EQ 0 &GOTO -NEXT 

x STYPE FILE &QI DATA A DOES NOT EXIST. USE ENTER TO EXIT. 

* &READ VAR S 

x & I F /& EQ / &EXIT 

X &GOTO -START 
-NEXT 

STATE RUNSAT JCL A 

& I F SRC EQ 0 &GOTO -COPY 

&TYPE TEMPLATE FILE RUNSAT JCL A DOES NOT EXIST ; 

&TYPE SEE LAB MANAGER. 

&EXIT -2 

-COPY 

TERROR &EXIT -1 

COPYFILE RUNSAT JCL A TEMP2 JCL A ( REPL 
& = &CONCAT OF / YOURQ1 FI L ENAME / SQ1 / 

&STACK TOP 
&STACK CHANGE S * 1 
&STACK FILE 
XEDIT TEMP2 JCL A 
EXEC SUBMIT TEMP2 JCL A 
CONWAIT 
DESBUF 
SEXIT 



Figure C.7 RUNSAT EXEC File 



FILE: RUNSAT JCL 



A1 



//RXHU5H JOB (4555,9999), , XCHAM , ,CLASS=A,REGI0N=1024K 
//G02 EXEC PGM=5ATL ITP 

//STEPLIB DD DS N= MSS. PHOENICS. LOAD, DISP=SHR 
// L I 5T1 DD DSN =MSS . PHOENICS . L I ST 1 .DIR/ DI5P=(0LD) 

//LIST2 DD DSN = MSS .PH0ENICS.LIST2.DIR,DISP=(0LD) 

//LIST3 DD DSN=MSS. PHOENICS. LIST3. DIR, DISP=(OLD) 

//LIST4 DD DSN=MSS. PHOENICS. LIST4. DIR, DISP=(OLD) 

//Q 1 DD DSN=MSS. SA555 .PHOENICS. SRC CYOURQ1 FI LENAME) , DISP=(OLD) 
//PHHELP DD DSN=MSS. PHOENICS. PHHELP. DATA, DISP=(OLD) 

//SCRL DD DSN=MSS. PHOENICS. SCRL,DISP=( NEW, PASS), 

// DCB=CRECFM=F,LRECL=80,BLKSIZE=80), 

// SPACE= ( TRK , (12,3) ), UNI T=SYSDA 

//COPY DD DSN=MSS. PHOENICS. COPY, DISP=( NEW, PASS), 

// DCB = (RECFM=F,LRECL=80,BLKSIZE=6320), 

// SPACE* (TRK, (12,3) ), UNI T=SYS DA 

//FT09F001 DD SYSOUT=* 

//FT06F001 DD SYSOUT=X 
//FT1AF001 DD SYSOUT=* 

//G512 DD DSN=MSS. PHOENICS. G512, DISP=(OLD) 

//G513 DD DSN=MSS. PHOENICS. G5 1 3 , DI SP= ( OL D ) 

//G51A DU DSN=MSS. PHOENICS. G514, DISP=(OLD) 

//G515 DD DSN=MSS. PHOENICS. G515, DISP=(OLD) 

//G516 DD DSN=MSS. PHOENICS. G516,DISP=(OLD) 

//G518 DD DSN=MSS. PHOENICS. G513 , DISP=( OLD) 

//G526 DD DSN=MSS. PHOENICS. G526 , DISP= ( OL D ) 

//G520 DD DS N= MSS. PHOENICS. G520,DISP=(OLD) 

//G521 DD DSN=MSS. PHOENICS. G521 , DI SP=( OL D) 

//G522 DD DSN=MSS. PHOENICS. G522 , DI SP= ( OL D) 

//G523 DD DSN=MSS. PHOENICS. G523, DISP=( OLD) 

//G52A DD DSN=MSS. PHOENICS. G52A , DI SP= ( OL D ) 

//G527 DD DSN=MSS. PHOENICS. G527 , DI SP=( OL D ) 

//G528 DD DSN=MSS. PHOENICS. G528 , DISP=( OLD) 

//G529 DD DSN=MSS. PHOENICS. G529, DISP=(OLD) 

//G530 DD DSN=MSS. PHOENICS. G530, DISP=(OLD) 

//G531 DD DSN=MSS. PHOENICS. G531 , DISP= ( OLD ) 

//G735 DD DSN=MSS. PHOENICS. G7 35 , D I SP= ( OL D ) 

//G736 DD DSN=MSS. PHOENICS. G7 36 , DI SP= ( OL D ) 

//G737 DD DSN=MSS. PHOENICS. G7 37 , DI SP= ( OL D) 

// DF1 0 DD DS N = MS S.S4555. PHOENICS . DF10, 

// DI SP=SHR 

//DF1 2 DD D S N = M S S . SA 555 . PHOENICS . DF1 2 , DI SP = SHR 



Figure C.8. RUNSAT JCL File 
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FILE: RUNEAR EXEC A1 


&TRACE 




EXEC SUBMIT RUNEAR JCL A 


FILE: RUNEAR JCL A1 


//EXHUSH 


JOB ( 4555 , 9999), ' CHAM ’ , CL ASS=G, REGION=2048K 


//xMAIN 0RG=NPGVM1 . 4555P, LINES* ( 12) 


//G02 EXEC 


PGM= EARTHP 


//STEPLIB 


DD DSN=MSS .PHOENICS . LOAD, DISP=SHR 


// DSTL 


DD DUMMY 


// DFO 9 


DD DSN=MSS . S4555 .PHOENICS. DF09 , DISP=SHR 


//DFOl 


DD DSN=&&DF01, DISP=( NEW, PASS), 


// 


DCB=(RECFM = VBS, LRECL = 19000, BLKSIZE = 19006 ) , 


// 


SPACE=(TRK, (25,3)),UNIT=SYSDA 


// DF02 


DD DSN=aSDF02, DISP=( NEW, PASS) , 


// 


DCB=(RECFM=VBS,LRECL=19000,BLKSIZE=19006), 


// 


SPACE=(TRK, (25,3) ),UMIT=SYSDA 


// DF03 


DD DSN=88DF03, DISP=( NEW, PASS) , 


// 


DCB=(RECFM=VBS,LRECL=19000,BLKSIZE=19006), 


// 


SPACE=(TRK, (25, 3) ),UNIT=SYSDA 


//DFO A 


DD DSH=88DF09,DISP=( NEW, PASS ) , 


// 


DCB=(RECFM = VBS, LRECL = 19000, BL KS IZE = 1 9 0 06 ) , 


// 


SPACE=( TRK, (25,3)), UMIT=SYSDA 


// DFO 5 


DD DSN = 88DF05,DISP=( NEW, PASS ) , 


// 


DCB = (RECFM=VBS, LRECL =190 00, BLKSIZE=19 006), 


// 


SPACE= (TRK, (25,3) ), UNI T=SYSDA 


// DF07 


DD DSN = 88DF07 , DISP = ( NEW, PASS ) , 


// 


DCB=(RECFM=VBS, LRECL =19000, BLKSIZE=19 006), 


// 


SPACE=(TRK,(25,3)),UNIT=SYSDA 


//DF08 


DD DSN=88DF08,DISP=( NEW, PASS), 


// 


DCB= (RECFM = VBS, LRECL = 19000, BLKSIZE = 190 06 ), 


// 


SPACE=(TRK, (25, 3) ),UNIT=SYSDA 


//DF15 


DD DSN=88DF15, DISP=( NEW, PASS) , 


// 


DCB=(RECFM = VBS, LRECL = 19 0 00, BLKSIZE = 19 0 06 ), 


// 


SPACE= (TRK, (25,3) ),UIIIT = SYSDA 


// DF16 


DD DSN=MSS . SA555 .PHOENICS . DF16 , DISP=SHR 


// DF1 0 DD 


i DSN = MSS.S<4555 .PHOENICS . DF10 , DISP = ( OLD) 


//DF12 DD 


i DSN=MSS.S4555.PH0ENICS.DF12,DISP=(0LD) 


//FT14F001 


DD SYSOUT = X 


// FT 06 FO 0 1 


DD SYSOUT = X 


// 





Figure C.9 RUNEAR EXEC And JCL Files 
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Figure C.10. Test Facility Side View (Y-Z Plane) 
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1. Test Facility Top View (X-Z Plane) 
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Figure C.12. Vertical Y-Z Plane Zones 



99 



21X3 HZ 



91NI3Z 



sin nz 



frlNUZ 

SIN 1 1Z 
21NI1Z 



UNIIZ 



T1X31Z 



CO CM 

CO C\J O O O «— l 

I I (_• I— I I— « | 

CO CO 2: 2: 2: CO 

—I — 1 — I — I 1 1 

X X X X X X 

\ ww 






XI 



Figure C.13. Horizontal X-Z Plane Zones 
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X GROUP 3. X-DIRECTION GRID SPECIFICATION 

* NX = THE NUMBER OF CELLS IN THE X DIRECTION (INTEGER) 
X XUL AST = DISTANCE IN THE X DOMAIN (METERS) 

* 

NX = XNCBL1 + XNCMD1 + XNCMD2 + XNCMD3+XNCBL2+XNCBL 3 + XNCEX1 
XULAST=XLBL1+XLMID1+XLMID2+XLMID3+XLBL2+XLBL3+XLEXT1 



* 

XFRACC1 )=-XNCBLl; 
XFRACC 3)=XNCMD1; 
XFRAC(5)=XNCMD2; 
XFRAC(7)=XNCMD3; 
XFRAC(9)=XNCBL2; 
XFRACC 11 )=XNCBL3; 
XFRAC(13)=XNCEX1; 
* 



XFRACC 2 )=XLBL1/(XNCBL1*XULAST) 
XFRACC A) =XLMID1/(XNCMD1*XUL AST) 
XFRAC(6)=XLMID2/( XNCMD2*XUL AST) 
XFRACC 8 )=XLMID3/(XNCMD3*XUL AST) 
XFRACC 10) =XLBL2/ (XNCBL2XXUL AST ) 
XFRACC 12)=XLBL3/CXNCBL3*XULAST) 
XFRACC 1 A) =XLEXT1/CXNCEX1*XUL AST) 



^ GROUP 4. Y-DIRECTION GRID SPECIFICATION 
x NY = THE NUMBER OF CELLS IN THE Y DIRECTION (INTEGER) 
X YVL AST = DISTANCE IN THE Y DOMAIN (METERS) 



* 



NY=YNCMI D+YNCBL 1+YNCBL2+YNCEX1+YNCEX2+YNCBL3 
YVLAST=YLMID+YLBL1+YLBL2+YL EXT1+YL EXT2+YL BL 3 
* 



YFRACC1)=-YNCBL1; 
YFRAC(3)=YNCMID; 
YFRAC(5)=YNCBL2; 
YFRAC(7)=YNCEX1; 
YFRACC9) =YNCEX2 ; 
YFRACC11 )=YNCBL3; 
X 



YFRACC 2 )=YLBL1/CYNCBL1*YVL AST) 
YFRACC 4 )=YLMID/(YNCMID*YVL AST) 
YFRACC 6 ) = YL BL2/ ( YNCBL 2*YVLAST ) 
YFRACC 8 )=YLEXT1/CYNCEX1*YVL AST) 
YFRACC 10) =YLEXT2/C YNCEX2*YVL AST) 
YFRACC 12) =YLBL3/(YNCBL3*YVLAST) 



x GROUP 5. Z-DIRECT ION GRID SPECI FICAT I ON 
x NZ = THE NUMBER OF CELLS IN THE Z DIRECTION (INTEGER) 
x ZWL AST = DISTANCE IN THE Z DOMAIN (METERS) 



x 



NZ=ZNCEX1+ZNCIN1 

ZHlAST=ZLEXT1+ZL 

X 



ZFRACC 

ZFRACC 

ZFRACC 

ZFRACC 

ZFRACC 

ZFRACC 

ZFRACC 

ZFRACC 



2 ) =ZL EXTl/ C ZNCEX1 *ZWL AST ) 
4)=ZLINT1/(ZNCIN1*ZWLAST) 
6)=ZLINT2/(ZNCIN2*ZWLAST) 

3) =ZLINT3/(ZNCIN3*ZWLAST) 
10)=ZLINT4/(ZNCIN4*ZWLAST) 
12 ) =ZL I NT5/ C ZNCIN5*ZWL AST ) 
14)=ZLINT6/(ZNCIN6*ZWLAST) 
16 ) =ZL EXT2/ ( ZNCEX2*ZWL AST ) 



+ZNCIN2+ZNCIN3+ZNCIN4+ZNCIN5+ZNCIN6+ZNCEX2 
INT1+ZLINT2+ZLINT3+ZLINT4+ZLINT5+ZLINT6+ZLEXT2 

1)=-ZNCEX1; 



ZFRACC 
ZFRACC 3) = ZNCIN1; 
ZFRAC(5)=ZNCIN2, 
ZFRACC7 ) 5 ZN C I N 3 ; 
ZFRACC 9 ) =ZNCIN4 ; 
ZFRACC 11 )=ZNCIN5 
ZFRACC 13)=ZNCIN6 
ZFRACC15)=ZNCEX2 



Figure C.14. Q1 MEthod-Of-Pairs Coding Sequence 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



GROUP 3. X-DIRECTION GRID SPACING 
CARTES = T 
NX = 9 



XULAST 

METHOD 


= 9. 
OF 


. 900E+00 

PAIRS USED FOR GRID 


SETTING 




XFRAC C 


1) 


= 


-1.000E+00 


; XFRAC 


C 


2) 


S 


5.051E-02 


XFRAC C 


3) 


= 


1 . 000E+00 


; XFRAC 


C 


4) 


= 


5.051E-02 


XFRAC C 


5) 


= 


1 . 000E+00 


; XFRAC 


c 


6) 


= 


5 . 051E-02 


XFRAC C 


7) 


= 


1 . 000E+00 


; XFRAC 


c 


8) 


= 


5.051E-02 


XFRAC C 


9) 


= 


1 . 000E+00 


; XFRAC 


( 


10) 


= 


5. 051E-02 


XFRAC C 


11) 




3. 000E+00 


; XFRAC 


c 


12) 


= 


1 .818E-01 


XFRAC C 


13) 


= 


1 . 000E+00 


; XFRAC 


c 


14) 


= 


2. 020E-01 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



GROUP A. Y- 
NY 

YVLAST = 1. 

METHOD OF 
YFRAC C 1) 



YFRAC C 
YFRAC C 
YFRAC C 
YFRAC C 



3) 

5 ) 

7) 

9) 



YFRAC C 11) = 



DIRECTION GRID SPACING 
16 

700E+01 

PAIRS USED FOR GRID SETTING 
= -4.000E+00 ; YFRAC C 2) = 

= 2 . 00 QE+QO ; YFRAC C A) = 

= A . OOOE+0O ; YFRAC ( 6) = 

= 2.000E+00 ; YFRAC C 8) = 

= 1.000E+00 ; YFRAC C 10) = 

3.000E+00 ; YFRAC C 12) = 



2.9A1E-02 
2.941E-02 
2.9A1E-02 
1. 132E-01 
1.265E-01 
1.176E-01 






GROUP 5. Z-DIRECTION GRID SPACING 
PARAB = F 
NZ = 28 

ZWL AST = 6.A90E+01 



METHOD OF PAIRS USED FOR GRID SETTING 



ZFRAC 


( 


1) 


= 


-1 . 000E+00 


; ZFRAC 


C 


2) 


= 


3. 082E-02 


ZFRAC 


C 


3) 


= 


5 . 000E+00 


; ZFRAC 


c 


4) 


= 


7.242E-02 


ZFRAC 


c 


5) 


= 


1 . 000E+00 


; ZFRAC 


c 


6) 


= 


4.931E-02 


ZFRAC 


c 


7) 


= 


2 . 0 0 QE+0 0 


; ZFRAC 


c 


8) 


“ 


1.156E-02 


ZFRAC 


c 


9) 


= 


6 . 000E+00 


; ZFRAC 


c 


10) 


= 


1.5A1E-02 


ZFRAC 


c 


11) 


= 


5.000E+00 


; ZFRAC 


c 


12) 


= 


5. 023E-02 


ZFRAC 


( 


13) 


= 


6 .000E+00 


; ZFRAC 


c 


14) 


= 


2.157E-02 


ZFRAC 


( 


15) 


s 
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Figure C.15. EARTH Me thod-Of -Pairs Output 
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Figure C.17. System Domain Grid (Y-Z Plane) 
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Figure C.18. H vs Y-Axis at IY=1 to IY=10 
For IX=1 And IZ=16 
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Figure C.19. W vs Y-Axis At IY=1 To IY=10 
For IX=1 And IZ=16 
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Figure C.20. H vs X-Axis At IX=1 To IX=5 
For IY=5 And IZ=16 
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Figure C.21. W vs X-Axis At IX=1 To IX=5 
For IY=5 And IZ=16 
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Figure C.22. H vs Y-Axis At IY=1 To IY=10 
For IX=1 And IZ=20 
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Figure C.23. W vs Y-Axis At IY=1 To IY=I0 
For IX=1 And IZ=20 



110 



PATCHC PLOTS 


,PROFIL, 


1, 5, 5, 5 


, 20, 20, 


1, 1) 


PLOT ( PL0T8 


, HI , 3.000E+05, 1 . 000E+06 ) 






TABULATION OF ABSCISSA 


AND ORDINATES. . . 






X 


HI 








2 . 50 0E-01 


6 .376E+05 








7 .50 0 E-0 1 


6 . 165E+05 








1 .250E+00 


5 . 89 0E+05 








1 .750E+00 


5.679E+05 








2 . 250E+00 


5.577E+05 








VARIABLE 


HI 








MINVAL= 


3 . 000E+05 








MAXVAL = 


1 . 000E+06 








CELLAV= 


5.937E+05 








1.00 + + . 








+ 


0.90 + 








+ 


0.80 + 








+ 


0.70 + 








+ 


0.60 + 








+ 


0.50 H — 


H 






+ 


0 . 90 + 




H 


“ft" 




0.30 + 








+ 


0.20 + 








+ 


0.10 + 








+ 


0.00 +....+., 








+ 


0 .1 


.2 .3 


. A ,5 .6 .7 


.8 .9 1 . 


0 


THE ABSCISSA 


IS X 


. MI N= 2.50E-01 


MAX= 2.25E+00 







Figure C.24. H vs X-Axis At IX=1 TO IX=5 
For IY=5 And IZ=20 
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Figure C.25. W vs X-Axis At IX=1 To IX=5 
For IY=5 And IZ=20 



112 



PATCHCPL0T19 


,PR0FIL, 1, 1, 13, 13, 21, 26, 1, 1) 


PLOTC PL0T14 


, HI , 5 . 000E+05 , 6 . 000E+05) 


TABULATION OF 


ABSCISSA AND ORDINATES... 


z 


HI 


1 . 000E+00 


5 . 929E+05 


A.350E+00 


5 . 535E+05 


9 . 050E+00 


5 . 653E+05 


1 .375E+01 


5 . 895E+05 


1 . 345E+01 


5 . 905E+05 


2 . 315E+01 


5 .631E+05 


VARIABLE 


HI 


MINVAL= 


5 . 0 00E+05 


MAXVAL = 


6 .000E+05 


CEL L AV = 


5.683E+05 


1.00 + + . . 




0.90 + 


+ 


0.80 + 




0.70 + 


M — " ~ 


0.60 + 




0.50+ ' 


iT^ - + 


0 .A0 hr - 


+ 


0.30 + 


+ 


0.20 + 


+ 


0.10 + 


+ 


0.00 




0 . 1 


.2 .3 .9 .5 .6 .7 .8 .9 1.0 


THE ABSCISSA 


IS Z . MIN= 1.00E+00 MAX= 2.31E+01 


****************** 



Figure C.26. H vs Z-Axis At IZ=21 To IZ=26 
For IX=1 And IY=13 
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Figure C.27. V vs Z-Axis At IZ=21 To IZ=26 
For IX=1 And IY=13 
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Figure C.28. H vs X-Axis At IX=1 To IX=5 
For IY=13 And IZ=23 
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Figure C.29. V vs X-Axis At IX=1 to IX=5 
For IY=13 And IZ=23 
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